Models for Low-Mass X-Ray Binaries in the Elliptical Galaxies 
NGC3379 and NGC4278: Comparison with Observations 



T. Fragos 1 , V. Kalogera 1 , K. Belczynski 2 , G. Fabbiano 3 , D.-W. Kim 3 , N. J. Brassington 3 , 
L. Angelini 4 , R. L. Davies 5 , J. S. Gallagher 6 , A. R. King 7 , S. Pellegrini 8 , G. Trinchieri 9 , S. 

E. Zepf 10 , A. Zezas 3 

tassosf ragos@northwestern . edu , vicky@northwestern . edu , kbelczyn@nmsu . edu , 
gfabbiano@cf a. harvard. edu, kim@cf a. harvard. edu, 
nbrassington@head. cf a. harvard. edu, angelini@davide .gsfc .nasa.gov, 
rld@astro.ox.ac.uk, jsg@astro.wisc.edu, ark@star.le.ac.uk, 
silvia. pellegrini@unibo . it , ginevra.trinchieri@brera. inaf . it , 
zepf @pa . msu . edu , azezas@cf a . harvard . edu 

ABSTRACT 

We present theoretical models for the formation and evolution of populations 
of low-mass X-ray binaries (LMXB) in the two elliptical galaxies NGC 3379 and 
NGC 4278. The models are calculated with the recently updated StarTrack code 



( IBelczynski et al.l 120071 ). assuming only a primordial galactic field LMXB popu- 



lation. StarTrack is an advanced population synthesis code that has been tested 
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and calibrated using detailed binary star calculations and incorporates all the im- 
portant physical processes of binary evolution. The simulations are targeted to 
modeling and understanding the origin of the X-ray luminosity functions (XLF) 
of point sources in these galaxies. For the first time we explore the population 
XLF down to lumi nosities of 3 x 10 36 ergs' 1 , as probed by the most recent ob- 



servational results (IKim et al.l l2006al ). We consider models for the formation 
and evolution of LMXBs in galactic fields with different CE efficiencies, stellar 
wind prescriptions, magnetic braking laws and initial mass functions. We iden- 
tify models that produce an XLF in excellent agreement with the observations 
both in shape and absolute normalization. We also find that the treatment of the 
outburst luminosity of transient systems remains a crucial factor for the deter- 
mination of the XLF since the modeled populations are dominated by transient 
X-ray systems. 

Subject headings: Stars: Binaries: Close, Stars: Evolution, X-rays: Binaries, 
Galaxies: Ellipticals 



INTRODUCTION 



A Low Mass X-ray binary (LMXB) is a Roche lobe overflowing, mass-transfering binary 
system with a compact object accretor, either a black hole (BH) or a neutron star (NS), and 
a low mass (> 1 M ) donor. Since the late 80's it has been suggested that LMXBs should 
exist in early type galaxies, E and SO, and that they might even dominate the X-ray emission 
( jTrinchieri and Fabbianolll985l ; lFabbianolll989l : IKim et al.l 119921 ). The stellar populations in 
these galaxies are typically old and homogeneous. Massive stars have already evolved to 
compact objects and LMXBs are probably the only sources with X-ray luminosities above 
10 36 ergs -1 . 

Uncontroversial detection of LMXBs in early type galaxies became possible only this las t 
decade with Chandra's increased angular resolution (jFabbiand 120061 ; ISarazin et al.l 120001 ) . 
The spectra of individual X-ray sources are consistent with those expected from LMXB 



models and the L MXBs observed in the Milky-way and M31 ([Humphrey and Buotd 12006 



Irwin et al.l 120031 ) . For many galaxies observed with Chandra, the XLFs have been derived 
and they can usually be fitted with a single or a broken power law. The d etections limit for 



these surveys is usually a few times 10 37 ergs 1 . IKim and Fabbianol (120041 ) derived XLFs for 



14 early type galaxies and they included completeness corrections. Each XLF is well fitted 
with a single power law with cumulative slope between -0.8 and -1.2. The composite XLF of 
these galaxies though is not consistent with a single power law. There is a prominent break 



- 3- 



at (5 ± 1.6) x 10 38 ergs 1 , close to the Eddington luminosity (L^d) of a helium accreting 
jreak might be hidden in the individual XLFs due to poor statistics (see 



NS-LMXB. This 



also Sarazin et al 



2000 



K^undu et al 



studies by IJeltema et al 



2002 



Jordan et al. 2004; Gilfanov 2004). Other recent 



(bo03h . ISivakoff et al.l (bo03h and I Jordan et all feooi ) suggested a 
break of the XLF at a higher luminosity (~ 10 39 erg s _1 ). The exact position and the nature 
of these breaks are still somewhat controversial, as the correct interpretation of the observed 
XLFs relies significantly on the proper completeness correction when looking at luminosities 
close to the detection limit and small number statistics at the high-end of the XLF. 



Recent Chandra observations (IKim et al.l l2006al ) have yielded the first low-luminosity 
XLFs of LMXBs for two typical old elliptical galaxies, NGC3379 and NGC4278. The detec- 
tion limit in these observations is ~ 3 x lO^ergs" 1 which is about an order of magnitude 
lower than in most previous surveys of early type galaxies. The observed XLFs of the two 
ellipticals extend only up to 6 x 10 38 ergs _1 and are well represented by a single power law 
with a slope (in a differential form) of 1.9 ± 0.1. 

When Chandra observations are compared with optical images from Hubble or other 
ground based telescopes, it is generally found that a significant fraction of the LMXBs are 
inside globular clusters (GC). On average 4-5% of the GCs in a given galaxy are associated 



with a LMXB (with L x >' 



10 37 erg s- 1 ) 



while the fraction of LMXBs located in GC, varies 



from 10% to 70% depending on the type of the galaxy and its GC specific frequency. At 
present the origin and the properties of these systems, both in GCs and the field, are not 
yet well understood. It has been noted that XLFs at high luminosities for each sub-group 



ered ( IFabbiano 



Kundu et al. 


2007; 


Jordan et al. 


2004; 


Sarazin et al. 


dFabbiano et al. 


2007 


) and 


Voss and Gilfanov 


(2007) 



independently found that the two XLFs (GCs and field sources) show significant differences 
at low luminisities (below 10 37 erg sec -1 ), pointing to a different LMXB formation mechanism 
in GCs. 

The natural question that arises is whether (i) all LMXBs were formed in GCs through 
dynamical interactions and some eventually escaped or some GCs dissolved in the field, 
or (m) field LMXBs were born in situ through binary evolution of primordial binaries. The 
formation rates associated with these two possibilities are not unders tood well en ough to give 
accurate predictions and use the relative numbers in the samples. IJuettl (120051 ) has shown 
that the observed relationship between the fraction of LMXBs found in GCs and the GC- 
specific frequency in early-type galaxies is co nsistent with the galactic field LMXB population 
being formed in situ. Similarly, llrwinl (120051 ) compared the summed X-ray luminosity of the 
LMXBs to the number of GCs in a galaxy; in the case of all LMXBs having formed exclusively 
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in GCs, the two should be directly proportional regardless of where the LMXBs currently 
reside. Instead he found that the proportionality includes an additive offset implying the 
existence of a LMXB population unrelated to GCs. 

In the past, semi-analytical theoretic al models have be e n intr oduced for the study of 
the LMXB population in early galaxies. IWhite and Ghoshl (119981 ) studied the connection 
between the star formation rates of normal galaxies, i.e. galaxies without an active nucleus, 
and the formation rate of LMXBs and millisecond pulsars, assuming that all LMXBs are 
formed from primordial binaries. Considering a time- dependent star formation rate, they 
showed that the general relativity timescales relevant to the evolution of primordial binaries 
to LMXBs and to millisecond pulsars, lead to a significant time delay of the peak in the 
form ation rate of these popula tions after the peak in the star formation rate. In a followup 
work lGhosh and White! (120011 ). using several updated star formation rate models, calculated 
the evolution of the X-ray luminosity of galaxies. They found that different star formation 
models lead to very different X-ray luminosity profiles, so the observed X-ray profiles can be 
used as probes of the star formation history. Finally they compared their models with the first 
Chandra deep imaging observa tions, to conclude t hat th ese first results were consistent with 
current star formation models. iPiro and Bildstenl (120021 ) argued that the majority of LMXBs 
in the field of elliptical galaxies have red giant donors feeding a thermally unstable disk and 
stay in this transient phase for at least 75% of their life. The very luminous X-ray sources 
(L x > 10 39 ergs _1 ) detected in Chandra surveys have been suggested to be X-ray binaries 
with highly super-Eddington mass inflow near the accreting component. In elliptical galaxies 
these objects have been suggested by King 20 02 to be micro- quasar-like, as these galaxies 
contain no high mass X-ray binaries ( King!Eo02 ). More recently Ivanova and Kalogera ( 2006 ) 
also argued that this bright end of the XLF is most likely dominated by transient LMXBs 
with BH accretors during outburst and it can be used to derive constraints on the BH 
mass function in LMXBs; they also showed that the standard assumption of a constant 
transient duty cycle (DC) across the whole population seems to be inconsistent with current 
observations. 



Semi-analytical p opulation s ynthesis (PS) models of LMXBs have also been constructed 
for late type galaxies. IWul (120011 ) created a simple birth-death model, in which the lifetimes 
of the binaries are inversely proportional to their X-ray luminosity, and calculated the XLFs 
of spiral galaxies. His models reproduce some features, such as the luminosity break in the 
observed XLFs of spiral galaxies. The position of this break depends on the star formation 
history of the galaxy, and he suggested that it can be used as a probe of the galaxy's merger 
history. 



The formation of LMXBs in GCs via dynamical interactions is less well studied, since 
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apart from the binary stellar eyo 



ution, one has to also take into account the complex cluster 



dynamics. iBildsten and Deloyei (120041 ) considered a semi-analytical model for accretion from 



degenerate donors onto NSs in ultracompact binaries and showed that binaries with orbital 
periods of 8-10 minutes and He or C/O white dwarf (WD) donors of 0.06-0.08 M naturally 
provide the primary slope (-0.8 for cumulative form) typically derived from XLFs of elliptical 
galaxies. Ultra-compact systems are predicted to form in the dense GC environment and 
have relatively short pe rsistent lifetim e s ( < 3 x 10 6 yr) but they form continuously through 
dynamical interactions. Ilvanova et al.l (120071 ) presented PS studies of compact binaries con- 
taining NSs in dense GCs. They used StarTrack as their PS modeling tool in addition to a 
simplified treatment for the dynamical interactions. Their models produced a mixed popu- 
lation of LMXBs with red giant and MS donors, and ultra-compact X-ray binaries; relative 
formation rates can be comparable but the different sub-populations have very different 
lifetimes. 

In this paper we investigate the plausibility of an important contribution to the XLFs 
of these two galaxies from a primordial galactic field LMXB population using advanced PS 
simulations. In §2 we describe briefly the physics included in our PS code and explain 
in detail the way we are constructing the modeled XLFs and the treatment of transient 
LMXBs. We discuss the results of our simulations in §3: the modeled XLFs from different 
models, a statistical comparison with the observed XLFs of the elliptical galaxies NGC3379 
and NGC4278, and an analysis of the dependence of the modeled XLF properties on the PS 
parameters. Finally in §4 we discuss the implication of our findings and the caveats of our 
methods. 



2. LMXB Population Models 



For the models presented in this study we focus on LMXBs formed in the galactic 
field as products of the evolution of isolated primordial binaries. The standard f ormation 
channel (IBhattacharya and van den Heuvellll99ll ; iTauris and van den Heuvelll2006l ) involves 
a primordial binary system with a large mass ratio; the more massive star evolves quickly 
to the giant branch and the system goes into a common envelope (CE) phase. During this 
phase, the less massive star, which is still dense and unevolved, orbits inside the envelope of 
the primary and is assumed to remain intact. The orbit of the system changes dramatically 
though, as orbital energy is lost due to friction between the the unevolved star and the 
envelope of the giant. Part of the lost orbital energy is used to expel the envelope of the 
giant star. The fraction of the lost orbital energy that is used to heat up the envelope 
of the giant star and finally expel it, defines the CE efficiency factor ckce- The CE phase 
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results in a binary system with an unevolved low mass main sequence (MS) star orbiting 
around the core of the massive star in a tighter orbit. The massive core soon reaches core 
collapse to form a compact object, either a NS or a BH and the binary orbit is altered due 
to mass loss and possible supernova kicks. If the binary does not get disrupted or merge in 
any of the stages described above, angular momentum loss mechanisms, such as magnetic 
braking, tides and gravitational wave radiation, will shrink further the orbit and the low- 
mass companion may evolve off the MS; The companion star eventually overflows its Roche 
lobe, transfering mass onto the compact object and initiating the system's X-ray phase. An 
alternative formation channel for NS-LMXBs is through accretion induced collapse of a WD 
accretor into a NS. These systems have generally very low X-ray luminosity, and do not 
affect the LMXB population in the luminosity range that we are interested in this paper. 



2.1. Synthesis Code: StarTrack 



We perform the simulations presented here with StarTrack (iBelczynski et al.l 12002 



20071 ). a advanced PS code that has been tested and calibrated using detailed mass trans- 
fer calculations and observations of binary populations, and incorporates all the important 
physical processes of binary evolution: (i) The evolution of single stars and non- interacting 
binary components , from ZAMS to remnant formation, is followed with analytic formulae 



( IHurley et al.l |2000| ). Various wind mass loss rates that vary with the stellar evolution- 
ary stage are incorporated into the code and their effect on stellar evolution is taken into 
account, (ii) Throughout the course of binary evolution, the changes in all the orbital prop- 
erties are tracked. A set of four differential equations is numerically integrated, describing 
the evolution of orbital separation, eccentricity and component spins, which depend on tidal 
interactions as well as angular momentum losses associated with magnetic braking, gravi- 
tational radiation and stellar wind mass losses, (in) All types of mass-transfer phases are 
calculated: stable driven by nuclear evolution or angular momentum loss and thermally or 
dynamically unstable. Any system entering the Roche lobe over-flow (RLOF) is assumed to 
become immediately circularized and synchronized. If dynamical instability is encountered 
the binary may enter a CE phase. For the modeling of this phase we use the standard 
energy balance prescription, (iv) The SN explosion is treated taking into account mass-loss 
as well as SN asymmetries (through natal kicks to NSs and BHs at birth). The distribution 
of the SN kick magnitudes is inferred fro m observed velocit ies of radio pulsars. For this 



project we use the distribution derived by lHobbs et al.l (120051 ) which is a single Maxwellian 



with a = 265kmsec~ 1 . It is however assumed that NS formation via electron capture or 
accretion induced collapse does not lead to SN kicks, (v) Finally, the X-ray luminosity of 
accreting binaries with NS and BH primaries (both for wind-fed and RLOF systems) is cal- 
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culated. For RLOF-fed systems there is a distinction made between persistent and transient 
(systems that undergo thermal disk instability), while wind-fed systems are always consid- 
ered as persistent X-ray sources. The mass-transfer is conservative up to the Eddington limit 
for persistent X-ray binarie s, while transients are allowed to have slightly super- Eddington 
luminosity (up to 3 x Lecm) (ITaam et al.lll997l ). In all cases we apply appropriate bolometric 
corrections (^boi) to convert the bolometric luminosity to the observed Chandra band. A 
much more detailed descripti on of all code elements, treatments of physical processes and 
implementation is provided in iBelczynski et al.l (120071 ). 



2.2. Model Parameters for NGC3379 and NGC4278 



In this study we focus on trying to understand the XLF characteristics of th e two ellip 



tical ga laxies NGC3379 and NGC4278, observed with Chandra and reported by iKim et al. 



( 12006al ). In the development of our models we incorporate our current knowledge about the 
characteristics of the stellar population in these galaxies (see Tabled]). The observationally 
determined parameters of the stellar populations, such as their age and metallicity , or their 
total stellar mass, are similar for NGC3379 and N GC4278. This allows u s to d evelop the 
same models in our simulations for both of them. iTerlevich and Forbes! (120021 ) estimated 



the ages and metallicities of 150 elliptical and late type spiral galaxies using published hight 
quality spectral line indices. For NGC3379 they are reporting an age of 9.3 Gyr and a 
metallicity of [Fe/H]=0.16, while for NGC4278 the corresponding values are 10.7 Gyr and 
[Fe/H]=0.14. The two galaxies have very similar optica l luminosity and assum ing the same 
light to mass ratio, they should also have similar masses. ICappellari et al.l (120061 ) used I-band 
observations from the Hubble Space Telescope to calculate the total stellar mass of the two 
galaxies and they found them to be 8.6 x 1O 1O M and 9.4 x 10 10 M o for N3379 and N4278 
respectively. The ratio of the integrated LMXB X-ray luminosity to the optical luminosity is 
4 times smaller for NGC33 79 which also has six times lower G C specific frequency compared 
to NGC4278 



sec 



Kim et al.ll2006al : lAshman and Zepj 119981 ). 



Table 1. Galaxy Properties 



Parameter 


NGC3379 


NGC4278 


References 


Distance (Mpc) 


10.57 


16.07 


Tonrv et al. (20011 


Age (Gyr) 


9.3 


10.7 


Terlevich and Forbes (20021 


Metallicity ([Fe/H]) 


0.16 


0.14 


Terlevich and Forbes (20021 


Mass (M Q ) 


8.6 x 10 10 


9.4 x 10 10 


Cappellari et al. (20061 


GC Specific Frequency 


1.2 


6.9 


Ashman and Zepf (19981 
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There is however a number of parameters in our models, for which we do not have any 
direct guidance from observations. We have no information about the star formation history 
of the two galaxies and thus we assume a ^-function like star formation episode at time t=0. 
Unknown are also the Initial mass function (IMF) and the distributions of orbital separation 
and eccentricity for the primordial binary systems. We adopt two different initial mass 
functions: Scalo/Kroupa and Salpeter, while fo r the distribution s of th e orbital properties 
we follow the standard assumptions described in iBelczynski et al.l (120021 ) . Other parameters 
that can affect the final LMXB population are the binary fraction of the host galaxy, the 
magnetic braking law adopted and the CE efficiency (o!ceJE The specific parameters we 
used to model the ellipticals NGC3379 and NGC4278 are listed in Table [2 



2.3. Models for the X-ray Luminosity Function 



In our models we keep track of all the binary properties, including the mass-transfer rate 
(M), as a function of time for populations of accreting NS and BH. We use the mass-transfer 
rates to identify the persistent and transient sources in our simulation results. Binaries 
for mass transfer rate higher than the critical rate M cr - lt for the thermal disk instability 
( IDubus et al.l Il999l ; iMenou et al.ll2002l ). are considered persistent sources and their X-ray 
luminosity (L x ) is calculated directly from the mass transfer rate as 



GM a M 
Vboie — , 



1 In our calculations, we combine q;ce and A into one CE parameter, where A is a measure of the central 
concentration of the donor. In the rest of the text, whenev er we mention the CE efficiency ckce, we practically 
refer to the product acE x A (see IBelczynski et al.l (|2007l ) for details.) 



Table 2. Model Parameters for NGC3379 and NGC4278 



Parameter 


Notation 


Value 


Star Formation 




8— function at t = 


Population Age 




9-10 Gyr 


Metallicity 


Z 


0.03 


Total Stellar Mass 


M* 


9 x 1O 1O M0 


Binary Fraction 


F bin 


50% 


IMF 




Scalo/Kroupa or Salpeter 


CE Efficiency 


ice 


20% - 100% 


Magnetic Braking 




Rappaport ct al. (19831 or Ivanova and Taam (2003) 
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where the radius of the accretor (R a ) is 10 km for a NS and 3 Schwarzschild radii for a BH, 
e gives a conversion efficiency of gravitational binding energy to radiation associated with 
accretion onto a NS (surface accretion e = 1.0) and onto a BH (disk accretion e = 0.5), and 
T^boi is a factor that converts the bolometric luminosity to the X-ray luminosity in the Chandra 
energy band ( 0.3 - 8 keV) . For RLOF accreting BH this conversion factor is estimated to 
0.8 jMiller et all bpplh while for RLOF accreting N S r^d = 0.55 Jpi Salvo et al. 



be n 


><- 


2002 


; ] 



Maccarone and Coppil 120031 ; iPortegies Zwart et al.l 12004 ) . The two correction factors 



e and 7? bol are applied in both persistent sources and transient sources in outburst. 

In the context of the thermal disk instability model, mass transferring binaries with 
M < M cr i t are considered transient sources, meaning that they spend most of their life in a 
quiescent state (T qu i scen t), at which they are too faint to be detectable, and they occasionally 
go into an outburst. The fraction of the time that these systems are in outburst (T outburst ) 
defines their DC: 

■j-jq ^outburst 

T -i- T . 

A outburst i ± quiescent 

Observations of Galactic LMXBs show that t ransient systems spend most of their life in the 
quiescent state, hinting at a DC below 20% (ITanaka and Shibazakilll996l ). 



The outburst luminosity and the DC of transient is not well understood and cannot be 
calculated from first principles. Instead we have to rely primarily on empirical constraints and 
simple theoretical arguments. In our analysis we consider a number of different treatments 
of these parameters for transients, which we describe in what follows. 

As a first approximation there has been suggested that transient LMXBs emit at 
their Eddington lumi n osity (L Edd ) when they are in outburst. In a different approach 
Portegies Zwart et al.l (120051 ) derived an empirical correlation between the outburst lumi- 
nosity of Milky Way transient LMXBs with BH accretors and their orbital period P: 



r/boie x mm 



2 x L Edd ,2 x L 



Edd 



P 

10h 



(2) 



We can generalize this relation to all transient LMXBs in galaxies other than our own, but 
we note that there has not been any observational work that shows that NS-LMXBs follow 
a similar trend. 

A more physical treatment is to assume that in the quiescent state the compact object 
does not accrete (or accretes an insignificant amount of mass) and matter from the donor is 
accumulated in the disk. In the outburst state all this matter is accreted onto the compact 
object emptying again the disk. Taking into account a lso that the X-ray luminosity probably 
cannot exceed L^dd by more than a factor of 2 (cf. iTaam et al.l 119971 ). we end up with a 
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definition of the outburst luminosity as: 



r] ho i€ x mm ( 3 x L E dd, 5 x 




(3) 



In the equation above, DC is unknown. iDobrotka et al.l (120061 ) studied accretion disk 
models for cataclysmic variables that are thought to experience the same thermal disk insta- 
bility (dwarf novae). They found a correlation between the DC of the system and the rate at 
which the donor star is losing mass M d . The exact relation of these to quantities depends on 
the values of the disk's viscosity parameters, but the general behavior can be approximated 
by: 



DC 



M, 



(4) 



crit 



Plugging eq.pj into eq.([3]) we eliminate the DC dependence and get an expression for the 
outburst luminosity of a transient system that depends only quantities which are directly 
calculated in our population modeling: 



L x = r?boie x min 3 x L 



Edd, 



R,M d 



(5) 



The accretion disk models by IDobrotka et al.l (120061 ) assume accretion onto a compact 
object with a hard surface and it is not obvious that the same results will apply for accre- 
tion onto a BH. In order to take into account all the available information, empirical and 
theoretical, about LMXB transient behavior, we treat BH and NS-LMXBs differently and 
define the outburst luminosity as: 



%oie x 



min [2 x L Edd , 3 x L E dd (1^)] , for BH acc. 



mm (2xL--. ™^ 



^Edd, 



RaM crit 



(6) 



for NS acc. 



We note that for BH-LMXBs, we adopt a single DC value for simplicity and lack of 
other information, although we have no clear physical reason to believe that all BH systems 
have the same DC. It is believed that the DC of BH systems is smaller than that of NS 
systems and on the order of 5% (ITanaka and Shibazakil Il996l ). We found that in all our 
models NS accretors greatly outnumber BH accretors. BH systems only have an important 



contrib ution at high luminosities, where the error bars in the observed XLFs (see iKim et al. 



2006al ) are too large to give us any tight constraints for our models. As we will see in section 
3.1, this treatment eq.([6]) gives us the best agreement with observations. 
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To construct the XLF we consider a snapshot of the whole population at the time we are 
interested in and we identify the LMXBs as transient or persistent. If a system is transient 
we decide whether it is in outburst or in quiescence according to its DC and either assign 
an outburst luminosity or discard the system as quiescent and hence too faint to contribute 
to the XLF. We then construct the XLF by calculating the cumulative X-ray luminosity 
distribution of the sources that are detectable (persistent and transient in outburst). We 
note that the age of the elliptical galaxies NGC3379 and NGC4278 and hence their LMXB 
population is known only to within ~ 1 Gyr. Consequently we cannot just choose a unique 
snapshot of the population. Instead we construct the XLF by considering the time window 
of 9 - 10 Gyr and considering time slices separated by 1 Myr. We construct the XLF at each 
of these timeslices and we take the average to represent the XLF that corresponds to the 
time window of interest. Doing so, we also improve the statistics of our model sample. 

It is computationally impossible to evolve enough binaries that would correspond to the 
total initial number of binaries in an elliptical galaxy (~ 10 9 binaries). For each model we 
evolve 10 6 binaries which takes about two months of CPU time on a modern processor. We 
then normalize to the total mass of the galaxy in question, taking into account the initial 
binary fraction and the initial mass function. 



3. Results 



3.1. Exploring the parameter space 



One of the implicit weaknesses of PS models is the large number of free parameters 
that one can vary and fine tune in order to get the desirable result. There are physical 
processes involved in the evolution of a binary system, such as stellar winds and magnetic 
braking, which are not fully understood. In this case various prescriptions are typically 
used to model them. Fortunately: (i) the result of interest to us in this study (XLF) is 
not sensitive to all of these model parameters; (ii) we can use some empirical knowledge 
from observations to constrain these parameters. We study a total of 336 models. In Table 
[3] we list 28 combinations of PS input parameters we study (CE efficiency, initial mass 
function, wind strength) and in Table H] we list the 6 different prescriptions we use for the 
determination of the DC and outburst luminosity of transient LMXBs. For each of the 
parameter comb ination from Table s [3] an d HI we try two prescr i ptions for the magnetic 



braking law, by llvanova and Taaml (120031 ) and iRappaport et al.l (119831). Parameter s not 



mentioned here are set as in the standard model considered in iBelczynski et al.l (120071 ). We 
name our models using a combination of a number from Table [3] which denotes the PS 
parameters of the models, a letter from Table H] which denotes the prescription we use for 
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Table 3. Population Synthesis Models: For each of the models listed below we app l ied tw o 

magnetic braking law prescriptions. In the rest of the text the llvanova and Taaml (120031 ) 
prescription w ill be denoted with the "IT" superscript at the end of the model name, while 
the iRappaport et al.l (119831 ) prescription with the "RVJ" superscript. 



Model a C E IMF r) wind 



1 


0.2 


Salpeter 


0.25 


2 


0.2 


Scalo/Kroupa 


0.25 


3 


0.2 


Salpeter 


1.0 


1 


0.2 


Scalo/Kroupa 


1.0 


5 


0.3 


Salpeter 


0.25 


6 


0.3 


Scalo/Kroupa 


0.25 


7 


0.3 


Salpeter 


1.0 


8 


0.3 


Scalo/Kroupa 


1.0 


9 


0,1 


Salpeter 


0.25 


10 


0,1 


Scalo/Kroupa 


0.25 


11 


0.1 


Salpeter 


1.0 


12 


0,1 


Scalo/Kroupa 


1.0 


13 


0.5 


Salpeter 


0.25 


11 


0,5 


Scalo/Kroupa 


0.25 


15 


0,5 


Salpeter 


1.0 


1(5 


0,5 


Scalo/Kroupa 


1.0 


17 


0.6 


Salpeter 


0.25 


18 


0.6 


Scalo/Kroupa 


0.25 


19 


0.6 


Salpeter 


1.0 


20 


0.6 


Scalo/Kroupa 


1.0 


21 


0.7 


Salpeter 


0.25 


22 


0.7 


Scalo/Kroupa 


0.25 


23 


0.7 


Salpeter 


1.0 


21 


0.7 


Scalo/Kroupa 


1.0 


25 


1.0 


Salpeter 


0.25 


26 


1.0 


Scalo/Kroupa 


0.25 


27 


1.0 


Salpeter 


1.0 


28 


1.0 


Scalo/Kroupa 


1.0 



Table 4. Treatment of Transient LMXBs: We tried 6 different prescriptions for the 
determination of the DC and outburst luminosity of transient LMXBs. 



Model Name 


ix,NS 


ix.BH 


DCns 


DCbh 


A 


eq.© 


eq.© 


eq.© 


5% 


13 


eq.© 


eq.© 


1% 


1% 


C 


eq.© 


eq.© 


7% 


7% 


D 


eq.© 


eq.© 


15% 


15% 


E 


^Edd 


^Edd 


10% 


10% 


F 


eq.© 


eq.© 


10% 


10% 



13 



the treatment of transient LMXBs and a superscript which denotes the magnetic braking 
law used. 

The analysis of several population simulations enables us to identify general behaviors 
of how the LMXB population is affected by changes to different parameters. In Figure [T] we 
show the modeled XLFs for 36 selected models, separated in 6 panels. In each panels all the 
parameters that determine the formation and evolution of LMXBs are kept constant, and 
we just change the modeling of transient systems (their DC and outburst luminosity). The 
observed XLFs of the two galaxies, NGC3379 and NGC4278, are also plotted for comparison. 

We find that the stellar wind strength alters significantly only the BH-LMXB population 
and thus only the high luminosity region of the modeled XLF. Stellar winds are important 
mainly in the evolution of massive stars which, depending on the wind prescription assumed, 
may lose a significant part of their envelope. Weaker stellar winds result in more massive 
pre-supernova cores and thus to the formation of more BHs. On the other hand, the ini- 
tial mass function affects the LMXB population globally. A flatter (Salpeter) initial mass 
function increases the overall number of sources with luminosities above 10 37 ergs _1 but 
still gives results consistent with the observations. The modeled XLF is more sensitive to 
the CE efficiency «ce- Smaller ckce values lead to more mergers among LMXB progeni- 
tors and therefore decrease the overall formation rate of LMXBs. If the CE efficiency is 
set to a very small value (oq e < 0.2), it becomes very difficult to form any LMXB at all 



(IKalogera and Webbink! Il998l ). while a high CE efficiency «ce > 0.6, clearly overproduce 
LMXBs. It is also clear that the adopted magnetic braking law greatly affects the charac- 
teristics of the model LMXB population. The XLFs of all 336 models can be seen in the 
online material supplemental to this paper. 

The ratio of persistent to transient LMXBs produced from our models, including systems 
both in outburst and in quiescent, is on the order of 1:20. So even if transient systems have 
a small DC (DC < 10%), their contribution to t he shape of the XL F is important an d som e 



times they may even dominate the population iPiro and Bildstenl (as discussed by 120021 ) . 
From this simple argument, we can understand that the shape of the modeled XLFs should 
be very sensitive to the DC and the outburst luminosity of transient systems. Figure [T] 
confirms this last statement, as we see that the larger variation in the shape of the modeled 
XLFs comes from different prescriptions in the treatment of transient systems and not from 
changing the binary evolution parameters. 
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Fig. 1. — Model XLFs for 36 selected LMXB population. In each panel all the PS parameters 
are kept constant, and only the modeling of transient systems changes. For comparison the 
observed XLFs of NGC3379 (dark red) and of NGC4278 (orange) are drawn. The modeling 
of the transient systems and their outburst characteristics can be more important than the 
usual PS parameters; Commonly used assumptions such as assigning the outburst luminosity 
of the transient LMXBs to be equal to Z>Edd (transient treatment E) lead to XLFs clearly 
inconsistent with the observations. 
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Table 5. Statistical comparison of the models shown in Figure (CQ) to the observed XLFs 
of NGC3379 and NGC4278. We used four different statistical tests for each case 
(Kolmogorov-Smirnov test, Kuiper test, Wilcoxon rank-sum test and x 2 goodness of fit 
test) and we are listing the probabilities, derived from each test, that a modeled XLF is 
consistent with the observed ones. For the statistical comparison of the complete list of 
models please see the electronic supplemental material. 



Model 




NGC 3379 






NGC 4278 






Pks 


-^Kuipcr 




Pks 


P K S 


^Kuiper 




Pns 


6A IT 


0.037 


0.065 


0.006 


0.022 


0.275 


0.020 


0.202 


0.046 


6B IT 


2.18e-04 


0.002 


o.ooe-Uo 


1.53c-04 


0.002 


0.007 


0.012 


3.92e-04 


6C IT 


3.89e-08 


1.10e-07 


7.75e-07 


1.10e-06 


4.89e-09 


*7 nT. nn 

7.^7e-U9 


5.96e-08 


1.04e-09 


6D IT 


1.61c-13 


1.68e-13 


U.UUc+UU 


5.48e-ll 


4.67e-15 


1. 12e-15 


n nr\„ 1 nn 

U.UUe+UU 


1.33e-15 


6E IT 


1.40e-24 


7.98e-25 


U.UUe+UU 


0.00e+00 


5.03e-34 


1.31c-35 


n nn „ 1 nn 

U.OOc+OU 


0.00e+00 


6F IT 


0.001 


l.bUe-Uo 


0.011 


0.00e+00 


6.60e-13 


1.97e-16 


n nn „ 1 nn 

U.UUc+UU 


7.46e-ll 


14A rvj 


2.77e-07 


D.59e-U7 


n nn „ i nn 

U.UUe+UU 


4.96e-10 


3.34e-04 


1.30c-04 


0.003 


6.39e-04 


14B rvj 


1.34e-10 


Q QQ^ i n 
O.OOC-IU 


n nrin i nn 
U.UUC+UU 


2.93c-10 


2.60e-08 


o.ooe-Uo 


0. loe-Uu 


4.68e-07 


14C RVJ 


4.94e-19 


1.83e-19 


n nn „ i nn 

U.OUc+OU 


0.00e+00 


3.61o-18 


2.11e-19 


n nn „ 1 nn 

U.OOc+UU 


0.00e+00 


14D RVJ 


6.87c-10 


c; c;o^ i n 
o.oze-lU 


U.UUe-f-UU 


0.00e+00 


6.96e-06 


l.zoe-Uo 


4.uoe-U4 


1.03c-10 


14E HVJ 


6.72c-27 


6.110-28 


0.00e+00 


0.00e+00 


5.05e-37 


2.14e-35 


0.00e+00 


0.00e+00 


14F RVJ 


4.38c-10 


3.35e-10 


2.38e-07 


0.00e+00 


3.17e-24 


3.97e-29 


0.00e+00 


0.00e+00 


1 A A IT 

14A 


0.804 


0.553 


0.351 


0.528 


0.661 


0.654 


0.302 


0.438 


14B IT 


0.780 


0.785 


0.291 


0.536 


0.890 


0.779 


0.462 


0.449 


14C IT 


0.049 


0.045 


0.171 


0.018 


3.24e-05 


4.32e-05 


3.29e-04 


3.36e-04 


14D IT 


1.70e-04 


3.22e-04 


0.001 


0.002 


1.14e-08 


5.67e-09 


0.00e+00 


2.76e-08 


14E IT 


3.82e-28 


3.38e-25 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


0.00e+00 


1 4F it 


1.95e-06 


5.16e-06 


0.001 


0.00e+00 


2.21e-08 


3.94e-09 


1.31e-06 


4.57e-07 


15A IT 


0.193 


0.303 


0.199 


0.065 


0.001 


0.004 


0.002 


2.36e-04 


15B IT 


0.871 


0.647 


0.366 


0.444 


0.837 


0.509 


0.483 


0.066 


15C IT 


0.298 


0.151 


0.454 


0.035 


2.92c-04 


5.57e-04 


9.97e-04 


5.96e-05 


15D IT 


0.052 


0.042 


0.194 


0.010 


2.04e-07 


2.11c-07 


2.86e-06 


2.87e-07 


15E IT 


1.72e-26 


1.51e-23 


O.OOc+00 


0.00e+00 


1.96e-44 


2.61e-41 


0.00e+00 


0.00e+00 


15F IT 


4.27e-06 


1.28e-06 


0.005 


l.llc-16 


4.95e-09 


2.47e-08 


1.19e-07 


1.67e-07 


16A IT 


0.832 


0.778 


0.294 


0.397 


0.498 


0.121 


0.377 


0.022 


16B IT 


0.585 


0.629 


0.204 


0.471 


0.276 


0.129 


0.328 


0.007 


16C IT 


0.009 


0.018 


0.044 


0.004 


5.22o-07 


3.61e-07 


3.91e-05 


6.95e-07 


16D IT 


1.91c-04 


3.26e-04 


0.001 


0.003 


1.54e-09 


1.16o-10 


0.00e+00 


8.82e-08 


16E IT 


8.10e-26 


5.55e-23 


0.00e+00 


0.00e+00 


6.45e-44 


3.10e-41 


0.00e+00 


0.00e+00 


16F IT 


2.78e-05 


9.86e-06 


0.008 


4.44e-16 


1.43e-10 


1.57c-13 


1.79e-07 


3.03e-09 


22A IT 


0.184 


0.566 


0.095 


0.606 


0.045 


0.103 


0.035 


0.358 


22B IT 


0.887 


0.921 


0.331 


0.793 


0.845 


0.575 


0.442 


0.426 


22C IT 


0.424 


0.122 


0.408 


0.046 


0.003 


7.28e-04 


0.025 


0.019 


22D IT 


0.041 


0.020 


0.160 


0.090 


3.17o-05 


9.40e-06 


2.01e-04 


6.87e-04 


22E IT 


6.17e-25 


2.35e-22 


O.OOc+00 


0.00e+00 


3.48e-43 


2.40e-40 


0.00e+00 


0.00e+00 


22F IT 


2.66e-04 


6.07e-05 


0.005 


2.66e-ll 


1.46e-09 


1.13e-08 


0.00e+00 


4.51e-06 
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We compare quantitatively our modeled XLFs to the observed ones from the two ellipti- 
cal galaxies NGC3379 and NNGC4278, using statistical tests that are not making an a priori 
assumption about the functional form of the observed XLF, and we report the probability 
that the two data sets are drawn from the same distribution. In order to make a fair com- 
parison between models and observations, we include in our statistical analysis only systems 
with luminosity above the completeness limit for the specific observations ( > 10 37 ergs~ 1 
for NGC3379 and > 3 x lCFergs" 1 for NGC4278). We have used three hypothesis testing 
statistical methods, the widely used Kolmogorov-Smirnov (KS), the Kuiper test which is a 
variation of the KS, equally sensitive in the whole range of the XLF and not only in the 
mid- range and the Wi lcoxon rank-sum test (RS). We have also used the x 2 goodness of fit 



test (IPress et al.l 120071 ). We note here that all four statistical tests are comparing only the 



shape of the XLFs and not their absolute normalization. 

In Table |5] we list the P values (probability that the two data sets are drawn from 
the same distribution) for all four tests, comparing the models shown in Figure [1] with 
each of the two observed XLFs from NGC3379 and NGC4278. For the statistical analysis 
of the complete list of models examined in this work, please see the online supplemental 
material. We use this analysis to draw general conclusions regarding the behavior of our 
models. In all cases, when we assume the outburst luminosity of the transient LMXBs to 
be equal to LEdd or use eq.© (transient treatment E and F, see Table 0J, the P values 
are extremely low and we can confidently say that this treatment of transient systems is 
highly unlikely to be correct. Assigning the same constant DC to all transients and then 
calculating the outburst luminosity using eq.®, sometimes leads to results consistent with 
the observations depending on the rest of the parameters used in the models. When this 
constant DC is as low as 1% (transient treatment B), the modeled XLFs resemble remarkably 
the observed ones. This happens because such a low DC practically removes the contribution 
of transient LMXBs to the total XLF. However, a DC of 1% is unrealistically low based on 
both observational evidence and theoretical predictions. Finally, the transient treatment A 
consistently provides the best agreement with observations, and in addition it is the one 
most physically motivated. One of our primary conclusions from this analysis is that careful 
modeling of the transient systems and their outburst characteristics is very important, more 
important than the usual PS parameters. 

Studying the complete list of models (see online supplemental material) becomes evi- 
dent that the magnetic braking law applied during the evolution of t h e LM XBs drastically 
changes the resulted population. Models where the iRappaport et al. I Jl983h pr escrip tion is 



used, produce XLFs inconsistent with the observations. The IRappaport et al.l (119831 ) brak 



ing law prescription pred icts much stronger angular momentum losses, compared to the 



Ivanova and Taaml (120031 ) one. Very strong angular momentum loss due to magnetic break- 
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ing lead to population of LMXBs where only wide binaries have avoided a possible merger. 
Our findings in favor of a mil der magnetic braking law (jlvanova and Taaml 120031 ) are in 
agreement with earlier work by iKalogera et all (119981 . and references there in) . 



Our statistical analysis tests only the shape of the XLF. However, within the uncer- 
tainties of the total mass of the two galaxies most of models we consider turn out to have 
the right normalization, and therefore this is not a strongly discriminating constraint. We 
note also that the PS model normalization can be very sensitive to certain binary evolu- 
tion parameters, such as the distribution of the mass ratios between the initial masses of 
the two binary components. Without a proper mult i- dimensional coverage of the parameter 
space and given the uncertainties in the total galaxy mass, we cannot use the total number 
of LMXBs observed as a formal constraint. Instead, we remain satisfied that most of the 
models we consider give us a total number of LMXBs, in the observed luminosity regime, 
with a facto r of 3 from the observe d sample. This factor is comparable to the galaxy mass 
uncertainty (ICappellari et al.ll2006l ) . 

In the discussion that follows we are using model 14A IT as our standard model that 
gives an XLF closest both in shape and normalization to the observed ones. 



3.2. Analyzing the LMXB population 

In order to investigate further the characteristics of our modeled LMXB population we 
classify the LMXBs based on their donor stellar type and also separate them in transient and 
persistent systems (Figure [2]). We also examine the orbital periods of each sub-population. 
This way we can infer to what degree each sub-population contributes to the total XLF. 
For our standard model 14A IT we find that systems with MS donors are the most numerous 
group, but their luminosity usually falls below the observational limit (~ 3 x 10 36 ergs _1 ). 
The XLF is dominated by transient and persistent systems with red giant donors, which 
despite being less numerous overall, they are more luminous. It is wor th pointing out that 



the po pulation of ultra-compact LMXBs with WD donors argued by iBildsten and Deloye 



(120041 ) to dominate the LMXBs in GCs, is also present in our models. Their semi- analytically 
derived XLF has, as expected, the same shape as the one found in our models (power-law 
with a slope of -0.8). We find in our models that the number of these systems formed in the 
galactic fields can be comparable to the number of LMXBs with red giant donors, but most 
of the time their contribution to the XLF is masked by the red giant donor systems. 



LMXBs with NS accretors outnumber those with a BH accretor by a factor of ~ 50. 
A small population of transient BH-LMXB has a significant contribution only to the high 
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luminosity end of the XLF, above 2 x 10 38 ergs -1 . These transient BH-LMXBs have or- 
bital periods from 1 to 10 hours and luminosity comparable to their Eddington luminosity 
according to eq.flH]) we adopted in this model. We also find a sub-population of persistent 
BH-LMXBs with MS donors, but their luminosity falls below the observational limit of the 
Chandra observations. 

We performed the same analysis for all of our models (plots similar to Figure [2] for all 
models can be found in the online supplemental material), and we identified some general 
trends on how the different PS parameters affect the LMXB population. The strongest effect 
comes from variation of the CE efficiency «ce- Low CE efficiencies (ckce < 0.3) result to the 
merger of all low period systems and only the widest binaries avoid a merger and survive 
the CE phase. Such wide binaries become wide orbit LMXBs with red giant donors. In 
models with higher CE efficiencies, binaries with progressively tighter orbits survive the CE 
phase, and we see LMXBs with shorter period and WD or MS donors appearing in our 
current population. Although these systems become more numerous for cvce > 0.4, LMXBs 
with red giant donors still have the most important contribution to the XLF for luminosities 
above 10 37 ergs -1 . 

The IMF affects the global picture of the LMXB population too. A flatter IMF 
(Salpeter) favors the formation of more massive stars. On the one hand this leads to the 
formation of BH-LMXBs and on the other hand favors the formation of generally more 
massive secondary stars which have evolved into WD by the time the system reaches the 
LMXB phase. We find that models with steep IMF (Scalo/Kroupa) have less LMXBs with 
WD donors than models with the same initial parameters but a flatter IMF (Salpeter). The 
stellar wind strength has similar effects on the LMXB population, as it affects the mass of 
the two components and thus their evolutionary state when they reach the LMXB phase. 
This effect is more prominent on the massive primary stars, where strong wind mass loss 
leads to less massive pre-supernova cores and thus fewer BH. 

Two of the PS parameters we varied in study (IMF and stellar wind strength) greatly 
affect the formation of BH-LMXB. Therefore some of our models (the ones with steep IMF 
and strong stellar winds) turn out to be highly inefficient in the formation of BH-LMXBs. 
However, a population of LMXBS with only NS accretors cannot form sources with lumi- 
nosities reaching up to 10 39 ergs -1 , unless we assume that highly super-Eddington accretion 
onto NS is possible. Consequently forming enough BH-LMXBs to populate the very high 
end of the XLF, is a discriminative criterion for our models. 
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Fig. 2. — Analysis of the LMXB population. We show the contribution of different sub- 
populations by separating the LMXBs into groups of systems with different donor (upper 
left panel) and accretor (lower left panel) stellar types. We find find that the mid-range of 
the XLF is dominated by NS-LMXB with red giant donors, while the high-end by BH LXBs 
with MS donors. Using the same separation into sub-populations we show the probability 
density function of the orbital period P of groups with different donor (upper right panel) 
and accretor (lower right panel) stellar types. 
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3.3. Time evolution of the LMXB population 

The age of the population is a very important factor in our simulations, as the charac- 
teristics of the binaries and their relative numbers and contributions vary significantly over 
time. In the right panel of Figure [3] we see the evolution of the modeled XLF versus time. 
The total number of sources as well as the number of luminous sources decreases steadily 
with time. At early times, the first 5 to 6 Gyr, the XLF has significant contributions from in- 
termediate mass X-ray binaries. The increased number of luminous sources at earlier times 
makes the shape of the XLF flatter. On the other extreme, in a 13 to 14 Gyr old stellar 
population, a population of luminous LMXBs is almost non-existent. The formation rate of 
X-ray binaries in general and of persistent sources in particular can be seen in right panel of 
Figure [3j The two curves follow closely each other, which means that according to our mod- 
els most the X-ray binaries go through a phase of persistent emission at some point of their 
evolution. Another characteristic that initially might seem counter-intuitive is that after a 
peak in the production of X-ray binaries in the first few gigayears, there is a decreasing but 
steady production of new persistent systems even after 10 or 11 Gyr. This suggest that there 
can be a very long delay between the formation of the compact object (at ~ 1 — 100 Myr) 
and the onset of th e mass-transfer phase . We s hould note that the derivation of the ages of 



the two galaxies by lTerlevich and Forbes! (120021 ) involves a systematic uncertainty of at least 



1 Gyr. This uncertainty in the age is an additional factor that can account for the difference 
in the number of observed sources between the two galaxies and of course some deviation 
between the observed and the modeled XLFs. 



3.4. The implication of GCs in the LMXB population: NGC3379 vs. 

NGC4278. 



The simulations we present in this initial study do not take into account formation 
channels involving dynamical interactions, believed to dominate the formation of LMXBs 
in GCs. However, we know from comparisons of optical and X-ray observations that 10- 
70% of LMXBs in early type galaxies are located in GCs, depending on the type of the 



galax y and the GC-specific frequency (IKundu et al.l 120071 ; I Jordan et al.l 12004 ; ISarazin et al. 
20031 ). It is thus clear that at least for some gala xies, the LMXB form ation in GC has 



an important contribution to the total population. iFabbiano et al.l (120071 ) studied the GC 
LMXB population of NGC3379 and found that only 15% of all LMXBs detected in this 
galaxy is coincident with GC. Furthermore, they argued that there are differences in the 
XLF of the two populations of field and GC LMXBs, but only at low luminosities (below 
10 37 ergs _1 ) where the incompleteness of the observed population does not allow reliable 
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direct comparison with theoretical models. Thus, it turns out that it is valid to compare our 
models to the total observed XLF (without subtracting the GC LMXBs) of NGC337, for 
luminosities above 10 37 erg s" 1 . The case of NGC4278 is different as its GC-specific frequency 
is 5 times larger than that of NGC3379, and the contr ibution of GC LMXBs is expected to 
be more significant there. An analysis similar to that of iFabbiano et al.l (120071 ) for NGC3379, 
will appear in the near future and it will enable us to separate the field from the GC LMXB 
population in our comparison. 

The realistic modeling of X-ray binaries formed i n GC requires the coupling of cluster dy- 
namics with binary evolution. iBildsten and Deloyd (120041 ) suggested that ultra-compact X- 
ray binaries dominate the cluster population and they are continuously replenished through 
dynamical interactions. These systems initially evolve through a bright persistent phase (for 
3Myr), then their mass transfer rate drops gradually as the WD donor mass decreases; the 
systems become transients and their contribution to the XLF is diminished. Our models 
for formation of LMXBs from primordial binaries produce a population of ultra-compact 
binaries too. Not surprisingly, since the properties of the binaries are very similar, the 
XLF of this sub-population has the exact same shape as the one produced by the model of 
Bildsten and Deloyd (120041 ): a power law shaped XLF (slope of -0.8), slightly flatter than the 
observed XLF of NGC3379 and NGC4278. In the case of the galactic field population though, 
there is no enhancement in the formation rate of these systems via dynamical interactions. 
So our models do not produce enough luminous persistent systems to dominate the total 
XLF. Their contribution, although significant, it is usually masked by LMXBs with red giant 
donors, which tend to have higher X-ray luminosities. More recently, llvanova et al.l (120071 ) 
studied the formation of compact binaries containing NSs in dense GCs. In their simulations 
they used Star Track as their PS model and adopted a simplified treatment of the dynamical 
events (see also llvanova and Kalogerall2006l ). They found that formation rate of dynamically 
formed LMXBs with red giant and MS donors is comparable to the ultra-compact X-ray bi- 
naries formed in GC, but the later ones have significant l y shor ter lifetimes. This last result 
casts doubts on the suggestion by IBildsten and Deloyd (120041 ) for an ultra-compact X-ray 
binary dominated population in GCs. 



Conclusions and Discussion 



The recent deep Chandra observations (IKim et al.ll2006al ) of the two typical old elliptical 
galaxies: NGC3379 and NGC4278, led to the first observed low-luminosity XLFs of LMXBs, 
with the detection limit (3 x 10 36 ergs _1 ) being about an order of magnitude lower than in 
most previous surveys. Motivated by this observational work, we developed PS simulations 
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of LMXBs appropriate for these two galaxies. We considered formation of LMXBs only 
through the evolution of primordial binaries in the galactic fields and examined the possible 
contribution to the overall LMXB pop ulation. For our modeling we used the updated PS 
code StarTrack ( IBelczynski et al.l 120071 ). 



Our main conclusions can be summarized as follows: 

We found that some of our models produce XLFs in very good agreement with the 
observations, based on both the XLF shape and absolute normalization. There is no unique 
combination of PS parameters and modeling of transient sources (DC and outburst lumi- 
nosity) that gives an XLF in agreement with the observation. We conclude that formation 
of LMXBs in the galactic field via evolution of primordial binaries can have a significant 
contribution to the total population of an elliptical galaxy, e specially the ones with low GC 
specific frequency such as NGC3379 OFabbiano et al.l 120071 ) . Nevertheless, we are able to 
exclude the majority of our models, as inconsistent with the observations. Note that widely 
used, simple assumptions such as that all transients source in the outburst state are emitting 
X-rays at LecM, ^ eac ^ ^° XLFs clearly inconsistent with the observed ones. Our results appear 
to be robust since we do not have to fine-tune our code parameters in order to get a model 
that resembles the observed population. 



As already suggested by iPiro and Bildstenl (120021 ). the LMXB population has a sig- 
nificant contribution by transient systems (thermal disk instability) and with reasonable 
outburst DCs they can even dominate the XLF. As a consequence the XLF shape is rather 
sensitive to the treatment of these transient systems. In Figure [1] we show that keeping 
the same PS parameters and changing only the modeling of transient sources leads to com- 
pletely different XLFs. We tried different methods of modeling the outburst characteristics 
of transient LMXBs and we found that: (i) When we assume the outburst luminosity of 
all transient LMXBs to be equal to LEdd (transient treatment E) or apply eq.([2]) (transient 
treatment F) - which was empirically derived for Galactic BH-LMXBs - to the whole pop- 
ulation, we get XLFs inconsistent with the observed ones regarding both their shape and 
the total number of sources predicted, (ii) A constant DC for all systems, although not 
physically motivated, can be sometimes a good first approximation, (in) We get the best 
agreement with obse rvations, when we con sider a variable DC for NS-LMXB based on the 
theoretical study of iDobrotka et al.l (120061 ) (see eqs. H]and[5]), while for the BH-LMXBs 
we use the empirical correlat ion between orbital period and outburst luminosity, derived by 
Portegies Zwart et al.l (120051 ) (see eq. |2J), and assuming a low constant DC (~ 5%). 



The LMXB sub-populations that mainly contribute to the model XLFs are NS-LMXBs 
with red giant donors and BH-LMXBs with MS donors. A population of persistent ultra 
compact LMXBs with WD donors is also present in our models and in some cases has an 
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important contribution too (see models 9A IT , 11A IT and 12A IT ). Of these sub-populations, 
the NS-LMXBs are the most dominant and primarily determine the XLF shape in the 
medium and low luminosity range (below 2 x 10 38 ergs _1 ), while the BH-LMXBs have a 
significant contribution to the high-end of the XLF. 

The normalization of the modeled XLFs is a less robust characteristic than its shape. We 
normalize the models so that the number of the primordial binaries we evolve, correspond 
to the known galaxy masses, given the IMF and the binary fraction. There are however 
uncertainties of the order of a few in the determination of the mass of the observed galaxies, 
due to uncertainties in their distance, the bolometric luminosity and the light-to-mass ratio. 
The majority of the models presented here produce the observed number of LMXBs to 
within a factor of 3, consistent with the galaxy mass uncertainties. Exceptions are models 
with transient treatment E or F (see Table H|) and models with high CE efficiencies (models 
21-28, see[3]) which greatly overproduce hight luminosity LMXBs. Furthermore small changes 
in the CE efficiency can change the total number of sources produced by a model by a factor 
of two, without changing significantly the shape of the XLF (compare for example models 
10A IT , 14A IT , 18A IT in the online supplemental material). In view of these uncertainties 
and our limited parameter space exploration for the models, we consider this normalization 
agreement satisfactory, but do not use it as an actual constraint on the models. 

We do not claim that the work presented in this Paper is a complete PS study of field 
LMXBs in elliptical galaxies. It is meant to be a first effort in interpreting the recent deep 
Chandra observations of the two elliptical galaxies NGC3379 and NGC4278. Throughout 
the Paper we identify the caveats of our analysis, and we intend to address them in our 
future work. It turns out that the realistic treatment of the outburst properties of transient 
LMXBs is crucial for the modeling of XLFs of extragalactic populations. The derivation 
of an empirical correlation between t he outburst luminosity and t he period of Galactic NS- 



LMXBs, similar to the one derived by lPortegies Zwart et al.l (120051 ) for Galactic BH-LMXBs 



will provide a better understanding of the differences in the transient behavior of these 
two classes of X-ray sources. Another simplifying assumption we made, is the constant 
outburst luminosity of transient LMXBs. The use of model lightcurves for the outburst 
phases of transient XLFs will possibly affect the shape of the modeled XLFs. We note that 
most of our models produce many systems with X-ray luminosity below the observational 
limit and down to 10 35 ergs~ 1 . The integrated diffuse X-ray emission from these galaxies 
can put a strong constraint on our models. The total luminosity of the observed diffuse 



emission will also include em ission from gas and stellar coronae ([Pellegrini and Fabbiano 



19941 ; iRevnivtsev et al.l 120071 ) and thus must be higher than the integrated luminosity of all 



the LMXBs in our models with luminosities below the observational detection limit. 
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Time (Gyr) L x (ergs/sec) 

Fig. 3. — Left panel: Formations rate of all x-ray binaries (blue line) and only of persistent 
x-ray binaries (red line) as a function of time. After a peak in the formation rate for the first 
few gigayears, there is a decreasing but steady production of new persistent systems even 
after 10 or 11 gigayears. Right panel: Evolution of the XLF with the age of the galaxy. 
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Statistical analysis and XLFs for the complete set of models (to appear as 

online supplemental material) 
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Table 6. Statistical Comparison of Models to the Observed XLF of NGC4278 



Model 




Pks 


NGC 3379 

^Kuiper ^X 2 


^RS 




Pks 


NGC 4278 

^Kuipcr ^X 2 


Prs 


1A RVJ 


2. 


85e- 


■11 


6 


19e-10 


5.96e-08 


1.33c-07 


9 


.52e-09 


1 


.74e-08 


3.76e-05 


6.50c-08 


1B RVJ 


7. 


41e- 


•12 


1 


.36e-10 


0.00e+00 


5.02c-08 


2 


,82e-09 


3 


.50c-08 


3.10e-06 


5.58c-08 


1C RVJ 


2 


51e- 


■17 


4 


.80e-17 


0.00e+00 


2.71c-14 


8. 


,12e-15 


6 


96e-15 


0.00e+00 


1.79c-ll 


1D RVJ 


3 


00c- 


-12 


4 


.18e-12 


0.00e+00 


2.79e-13 


2 


,73e-09 


4. 


99e-10 


O.OOe+OO 


2.39c-08 


1E RVJ 


9 


61e- 


■22 


7 


.91e-22 


0.00e+00 


O.OOc+00 


1 


46c-20 


1 


,78c-20 


0.00e+00 


O.OOc+00 


-j^pRVJ 


3 


30c- 


■20 


3 


.12c-20 


0.00e+00 


0.00e+00 


1 


,15c-19 


1 


,82c-22 


0.00e+00 


0.00e+00 


1A IT 


2. 


,77e- 


■10 


7 


57e-10 


0.00e+00 


7.93c-10 


2. 


,63e-04 


1 


21e-04 


5.53e-04 


1.34c-05 


1B IT 


1 


59e- 


■09 


1 


.10e-08 


0.00e+00 


5.23c-07 


7. 


,33e-05 


2 


lle-05 


3.50e-04 


6.18c-07 


1C IT 


2 


83e- 


■15 


4 


.29e-15 


0.00e+00 


2.89e-14 


1 


.08e-12 


6 


79e-13 


0.00e+00 


3.17c-13 


1D IT 


2 


OOe- 


■18 


1 


.14e-18 


O.OOc+00 


O.OOc+00 


6 


,23c-16 


1 


,71c-17 


O.OOc+00 


0.00e+00 


1E IT 


1 


37e- 


■19 


8 


.10c-20 


0.00e+00 


O.OOc+00 


5. 


llc-19 


6 


14e-20 


0.00e+00 


0.00e+00 


1F IT 


3 


03e- 


-12 


3 


.81e-12 


0.00e+00 


1.71e-12 


2 


,42e-15 


2 


00e-19 


0.00e+00 


5.00e-15 


2A RVJ 


1 


61e- 


■07 


1 


,68c- 06 


1.19e-07 


5.56c-06 


7 


,72e-05 


8. 


85e-05 


6.56e-04 


1.39e-04 


2B RVJ 


5 


81e- 


11 


6 


.54e-10 


0.00e+00 


1.86c-07 


4 


,83e-08 


3 


56e-07 


2.92c-06 


7.14e-07 


2C RVJ 


4. 


80e- 


■20 


4 


.32c-20 


O.OOc+00 


O.OOc+00 


5 


,02c- 17 


1 


.57e-17 


O.OOc+00 


2.22c-16 


2D RVJ 


1 


34e- 


■18 


6 


.52e-19 


O.OOc+00 


O.OOc+00 


2. 


,20c-14 


3 


33e-15 


0.00e+00 


4.44e-16 


2E RVJ 


2. 


55e- 


■22 


1 


.29c-22 


0.00e+00 


O.OOc+00 


9 


,22e-24 


3 


38e-23 


O.OOc+00 


0.00e+00 




2 


55e- 


■22 


1 


.29c-22 


0.00e+00 


O.OOc+00 


2 


12c-21 


4. 


96c- 24 


O.OOc+00 


O.OOc+00 


2A IT 


2 


14e- 


09 


9 


,96c-09 


0.00e+00 


1.27c-06 


8 


.22e-05 


3 


19e-05 


3.37e-04 


5.67c-06 


2B IT 


2. 


14e- 


■09 


2 


.08c-08 


0.00e+00 


8.51c-07 


7 


.23e-07 


4 


50c-06 


7.34c-05 


4.87e-08 


2C IT 


4. 


07e- 


■14 


8 


,46c-14 


0.00e+00 


6.49c-ll 


1 


,64c-13 


1 


43e-13 


O.OOc+00 


5.05c-13 


2D IT 


3 


32e- 


■18 


2 


.73e-18 


0.00e+00 


0.00e+00 


1 


,73c-16 


4 


23e-17 


O.OOc+00 


3.33c-16 


2E IT 


3 


,22e- 


■21 


1 


.78e-21 


0.00e+00 


0.00e+00 


2. 


22e-22 


6 


42e-23 


O.OOc+00 


O.OOc+00 


2F IT 


1 


,37e- 


■08 


2 


.64c- 08 


0.00e+00 


3.70c-ll 


4 


41c-12 


8. 


13e-15 


O.OOc+00 


2.83e-10 


3A RVJ 


2 


04e- 


-06 


3 


.09c-05 


8.17e-06 


4.20c-06 


2 


.69e-04 


3 


62e-04 


0.016 


3.45e-05 


3B RVJ 


1 


18e- 


■05 


7 


,48e-05 


1.08e-04 


1.46e-04 


7 


.87e-04 


3 


92e-05 


0.050 


1.87c-06 


3C RVJ 


4. 


08e- 


13 


2 


.72e-12 


0.00e+00 


1.09e-09 


2 


.02e-10 


8. 


26e-10 


O.OOc+00 


7.45e-10 


3D RVJ 


3 


,24e- 


■15 


6 


55e-15 


O.OOc+00 


5.99c-13 


1 


.llc-12 


1 


.01c-12 


0.00e+00 


1.19e-10 


3E RVJ 


5 


05e- 


■16 


1 


.56c-15 


O.OOc+00 


1.57c-12 


2 


,04c-13 


1 


.05c-12 


O.OOc+00 


4.40e-14 


3F rvj 


1 


03e- 


■14 


3 


.48c-14 


O.OOc+00 


3.69c-ll 


4 


42c-13 


2 


.10c-14 


0.00e+00 


1.39c-10 


3A IT 


3 


16e- 


■07 


1 


■87c-06 


5.96c-08 


6.86c-07 




0.002 




0.002 


0.003 


0.002 


3B IT 


3 


98e- 


■06 


6 


.67c-05 


5.60c-06 


7.51e-05 




0.001 


6 


.llc-04 


0.009 


1.12e-04 


3C IT 


4. 


,74e- 


■12 


1 


,78c- 11 


0.00e+00 


8.18c-09 


1 


,29c-10 


2. 


.80c-10 


O.OOc+00 


2.68e-10 


3D IT 


2. 


.16e- 


■16 


2 


,88c-16 


0.00e+00 


7.55e-15 


1 


99c-13 


2. 


.72c-14 


0.00e+00 


1.29c-ll 


3E IT 


1 


96e- 


■19 


1 


■76c-19 


O.OOc+00 


O.OOc+00 


2. 


63e-20 


2 


,80e-20 


0.00e+00 


0.00e+00 


3F IT 


5. 


13e- 


■07 


1 


.36e-06 


l.Olc-06 


3.29c-08 


7 


.66c-12 


2 


,21c-13 


O.OOc+00 


1.24c-09 


4A R VJ 


1 


47e- 


■10 


1 


.66e-09 


0.00e+00 


4.10c-08 


1 


47c-07 


1 


.18c-07 


1.64c-05 


3.95e-07 


4B R VJ 


1 


,47e- 


■11 


2 


49c-10 


O.OOc+00 


3.11c-08 


2 


.81c-08 


3 


16c-07 


7.51c-06 


2.03c-07 


4C RVJ 


4. 


,07c- 


-18 


6 


51c-18 


O.OOc+00 


l.lle-16 


5. 


,82c-15 


4 


.08c-15 


0.00e+00 


1.39c-13 


4D RVJ 


4. 


23e- 


■20 


2 


,47c-20 


O.OOc+00 


O.OOc+00 


3 


.58c-17 


6 


.00c-18 


O.OOc+00 


O.OOc+00 


4E RVJ 


6 


Ole- 


■21 


5 


.06e-21 


O.OOc+00 


0.00e+00 


5 


,08e-20 


1 


,12e-20 


O.OOc+00 


0.00e+00 


4F rv.i 


6 


Ole- 


■21 


5 


,06c-21 


O.OOc+00 


O.OOc+00 


4 


34c-19 


6 


35e-22 


O.OOc+00 


0.00e+00 


4A IT 


1 


70c- 


■10 


6 


78c-10 


O.OOc+00 


1.72c-09 


1 


91c-07 


2 


80e-08 


1. llc-04 


l.Olc-07 


4B IT 


5. 


39e- 


13 


3 


,17c-12 


O.OOc+00 


1.64e-09 


2 


66c-10 


9 


.41c-10 


2.24c-05 


9.39c-ll 


4C IT 


7. 


16e- 


■19 


6 


.41e-19 


0.00e+00 


l.lle-16 


1 


,75e-18 


4 


19c-19 


O.OOc+00 


O.OOc+00 
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Table 6 — Continued 



Model 




NGC 3379 






NGC 4278 






Pks 


-^Kuiper ^~X 2 


^RS 


■F>KS 


^Kuiper ^~X 2 


PrS 



4D IT 


2.93e-21 


4E IT 


1.12e-23 


4F IT 


2.58e-12 


5A rvj 


7.51e-09 


5B RVJ 


3.27c-09 


5C RVJ 


1.62e-16 


5D RVJ 


8.30e-17 


g E RVJ 


l.lOc-22 


gpRVJ 


2.17e-14 


5A IT 


0.002 


5B IT 


5.11e-04 


5C IT 


2.09c-04 


5D IT 


2.42e-05 


5E IT 


1.08e-22 


5F IT 


0.006 


6A rvj 


4.97e-13 


6B RVJ 


9.71c-14 


g C RVJ 


6.86c-19 


6D RVJ 


3.70e-18 


6E RVJ 


1.16e-23 


gpRVJ 


7.06c-16 


6A IT 


0.037 


6B IT 


2.18e-04 


6C IT 


3.89c-08 


6D IT 


1.61c-13 


6E IT 


1.40e-24 


6F IT 


0.001 


7A RVJ 


4.18e-15 


7B RVJ 


2.72e-15 


7C RVJ 


3.41e-22 


7D RVJ 


1.45e-23 


7E RVJ 


6.10e-25 


7pRVJ 


1.04e-16 


7A IT 


1.25c-04 


7B IT 


3.57c-06 


7C IT 


3.41e-06 


7D IT 


5.07c-07 


7E IT 


5.72e-24 


7F IT 


6.45e-04 


8A rvj 


8.95e-17 


g B RVJ 


2.30c-15 


g C RVJ 


1.48e-22 


g D RVJ 


1.27c- 24 


8E RVJ 


3.24e-25 


gpRVJ 


1.45e-18 



1.13e-21 
3.55e-24 
2.76c-12 
8.03c-08 
6.74e-08 
2.82c-16 
5.68c-17 
5.47c-23 
2.44e-14 
0.004 
0.003 
7.37c-04 
6.19c-05 
4.31e-20 
3.04e-06 
1.81e-12 
5.53e-13 
6.57e-19 
1.93c-18 
3.93c-24 
5.44e-16 
0.065 
0.002 
l.lOc-07 
1.68c-13 
7.98c-25 
1.60e-05 
2.33c-14 
1.44e-14 
2.16c-22 
3.82c-24 
1.61c-25 
6.62c-17 
9.19c-05 
3.03c-05 
1.35e-05 
1.28c-06 
4.54e-23 
3.79c-07 
3.67c-16 
1.47c-14 
8.35c-23 
2.89c-25 
8.42c-26 
7.66e-19 



0.00e+00 
O.OOc+00 
0.00e+00 
1.19c-07 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
0.00e+00 
7.23e-04 
2.38c-04 
1.59e-04 
3.16c-04 
O.OOc+00 

0.109 
0.00e+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 

0.006 
5.63e-05 
7.75c-07 
O.OOc+00 
O.OOc+00 

0.011 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
0.00e+00 
5.56c-05 
7.75c-07 
1.91c-06 
7.51e-06 
O.OOc+00 

0.041 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOe+00 



0.00e+00 
O.OOc+00 
O.OOc+00 
5.51e-06 
3.04e-06 
4.83c-14 
3.33c-16 
0.00e+00 
l.lle-16 
0.001 
0.004 
6.43c-05 
2.68c-06 
O.OOc+00 
7.62c-ll 
1.88c-09 
1.03c-09 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 

0.022 
1.53c-04 
1.10c-06 
5.48c-ll 
O.OOc+00 
O.OOc+00 
8.46e-ll 
5.54c-ll 
O.OOc+00 
O.OOc+00 
O.OOc+00 
O.OOc+00 
1.09c-04 
1.00c-05 
2.07c-07 
9.61e-06 
O.OOc+00 
O.OOc+00 
1.89c-12 
5.68c-ll 
O.OOc+OO 
O.OOc+00 
O.OOc+00 
0.00e+00 



1.70e-19 
3.75c-23 
1.04c-17 
4.97c-07 
7.86c-07 
4.51c-14 
2.33c-15 
6.92e-28 
3.63c-20 
0.152 
0.014 
1.64e-04 
1.49c-06 
3.12e-31 
9.75c-08 
8.67c-ll 
1.26c-10 
2.53c-16 
2.41c-12 
2.58e-22 
2.34c-20 
0.275 
0.002 
4.89c-09 
4.67c-15 
5.03e-34 
6.60e-13 
1.31c-ll 
1.33c-12 
4.41c-20 
7.48e-23 
2.00e-24 
5.64e-23 
0.092 
0.002 
6.38e-05 
1.63e-07 
2.58e-31 
4.96e-09 
3.91c-14 
1.17c-12 
4.33e-21 
3.56c-24 
2.08c-25 
1.18e-23 



1.27c-20 
1.80e-24 
1.35e-21 
1.96e-06 
1.06c-05 
3.30c-14 
7.93c-17 
2.63e-28 
3.15e-23 
0.004 
0.005 
1.22c-04 
8.33c-07 
1.82c-28 
4.92c-10 
6.62c-12 
4.37c-10 
8.14c-17 
5.46c-13 
1.56e-23 
9.01c-24 
0.020 
0.007 
7.27c-09 
1.12c-15 
1.31c-35 
1.97c-16 
1.37c-ll 
4.09c-12 
7.57e-21 
3.33e-25 
7.39c-26 
8.06e-27 
0.002 
0.001 
7.95e-05 
8.30e-08 
3.53c-31 
1.49c-ll 
2.41c-15 
4.46c-12 
5.77c-22 
l.Olc-25 
6.83e-27 
1.51e-28 



0.00e+00 
0.00e+00 
0.00e+00 
4.83e-06 
2.32e-06 
O.OOc+00 
0.00e+00 
O.OOc+00 
O.OOc+00 
0.156 
0.032 
0.005 
5.96c-07 
0.00e+00 
5.60c-06 
1.55e-06 
1.79c-06 
0.00e+00 
O.OOc+OO 
0.00e+00 
0.00e+00 
0.202 
0.012 
5.96c-08 
0.00e+00 
0.00e+00 
0.00e+00 
1.25e-06 
1.19c-07 
0.00e+00 
0.00e+00 
O.OOc+OO 
0.00e+00 
0.243 
0.011 
0.002 
5.96e-08 
0.00e+00 
1.73c-06 
5.96c-08 
4.77c-07 
O.OOc+OO 
O.OOc+00 
0.00e+00 
O.OOc+OO 



0.00e+00 
0.00e+00 
O.OOc+OO 
4.00c-06 
9.27c-07 
2.31c-12 
0.00e+00 
0.00e+00 
O.OOc+OO 

0.002 
6.97c-04 
5.43e-04 
5.43e-07 
O.OOe+00 
3.73e-06 
1.35e-09 
6.94e-10 
7.22c-15 
O.OOe+00 
O.OOc+OO 
O.OOc+00 

0.046 
3.92e-04 
1.04e-09 
1.33c-15 
O.OOe+00 
7.46e-ll 
2.48e-ll 
1.20c-ll 
O.OOc+OO 
O.OOe+00 
O.OOe+00 
O.OOe+00 

0.001 
8.81c-05 
1.20e-04 
3.18e-08 
O.OOe+00 
3.23c-07 
8.46e-14 
3.77c-13 
O.OOc+OO 
O.OOc+00 
O.OOe+00 
O.OOc+OO 
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Table 6 — Continued 



Model 




NGC 3379 






NGC 4278 






^KS 


Pa- ■ 

1 Kuipcr 


P •> 




PkS 


1 Kuiper 


P Q 




o A IT 

8 A 1 L 


8.67e-06 


2.95e-05 


1.43c-06 


2.97e-05 


0.032 


0.006 


0.014 


7.16e-04 


8B ±X 


2.10e-05 


1.72c-04 


2.28e-05 


0.001 


5.90c-04 


6.00c-05 


0.005 


1.37e-06 


8C 


4.26e-12 


1.00e-ll 


0.00e+00 


7.34e-09 


7.66e-12 


8.02c-12 


O.OOc+00 


7.25c-ll 


ot^TT 

8D 11 


5.40e-15 


4.67e-15 


O.OOc+00 


4.54c-14 


2.72e-14 


1.29c-15 


0.00e+00 


3.22e-15 


ot-'IT 

8E 1J - 


1.72e-25 


2.40e-24 


0.00e+00 


0.00e+00 


3.70e-35 


1.06c-34 


O.OOc+00 


0.00e+00 




0.002 


5.24e-07 


0.045 


4.44e-16 


8.17e-13 


2.73c-16 


O.OOc+00 


3.57e-ll 


n A RV I 
9A^ J 


8.54e-10 


3.37c-09 


4.77c-07 


5.22e-07 


4.45e-10 


5.07c-ll 


1.55e-06 


1.21e-09 


hdRVJ 
9B 


1.52e-10 


7.08c-10 


1.19e-07 


2.09e-07 


1.05c-10 


3.28e-10 


1.19c-07 


4.84c-10 


9C 


8.77e-17 


8.58c-17 


0.00e+00 


6.72e-14 


8.09c-16 


1.64c-16 


0.00e+00 


3.33e-16 




6.60c-13 


5.36c-13 


0.00e+00 


2.61e-13 


2.02c-12 


3.84c-13 


0.00e+00 


6.32e-12 


9E 


1.90e-24 


3.82e-25 


O.OOc+00 


0.00e+00 


8.02e-31 


1.44c-28 


0.00e+00 


0.00e+00 


9F 


2.49e-09 


2.94c-09 


1.79c-06 


2.22e-16 


3.33c-22 


5.60e-26 


O.OOc+00 


0.00e+00 


r> A IT 


0.552 


0.161 


0.327 


0.054 


0.006 


0.007 


0.025 


0.006 


9B 


0.729 


0.751 


0.410 


0.519 


0.103 


0.075 


0.046 


0.070 


n/"iTT 

9C 11 


0.225 


0.095 


0.121 


0.019 


0.006 


0.017 


0.029 


0.005 


nT\IT 


0.164 


0.010 


0.447 


7.73e-04 


2.08e-07 


1.27c-07 


9.30c-06 


2.25e-06 


nr->TT 
9E 11 


1.31c-26 


7.56e-24 


0.00e+00 


0.00e+00 


1.40e-44 


1.17c-41 


O.OOc+00 


O.OOc+00 




8.81c-07 


8.34c-07 


0.001 


7.77e-16 


4.60c-06 


1.52c-06 


2.84e-04 


2.60c-05 


1 n A RV 7 


1.31&-05 


6.61e-05 


3.52e-06 


3.62e-05 


1.04e-04 


2.12c-04 


0.001 


7.43c-04 


-i nD RV 7 


9.46e-07 


6.46c-06 


1.97c-06 


6.71e-06 


6.70e-06 


3.46e-05 


7.34e-04 


1.48e-05 


i nriRV I 
IOC 


6.67e-17 


6.68c-17 


O.OOc+00 


2.55c-15 


7.62e-17 


2.07c-17 


O.OOc+00 


9.99c-16 


i nnRV T 


2.51e-14 


1.68c- 14 


O.OOc+00 


0.00e+00 


1.04e-08 


3.52e-09 


O.OOc+00 


4.95e-12 


1 nnRV I 

10E 


2.88e-24 


5.86e-25 


O.OOc+00 


O.OOc+00 


6.47e-33 


9.84c-35 


0.00e+00 


O.OOc+00 


i nrRV 7 

10F^ J 


8.43e-08 


1.15e-07 


9.78c-06 


l.llc-16 


7.51e-21 


1.82e-24 


O.OOc+00 


0.00e+00 


1 n A TT 


0.073 


0.156 


0.072 


0.144 


0.003 


5.77c-04 


0.021 


0.004 




0.265 


0.480 


0.089 


0.650 


0.043 


0.014 


0.044 


0.004 


IOC 11 


9.08e-04 


3.06e-04 


0.035 


3.22c-04 


1.15e-10 


6.54c-ll 


0.00e+00 


9.64e-10 


i nnTT 


7.o7e-09 


8.75c-09 


5.01c-06 


2.09c-07 


5.24e-17 


5.08e-18 


O.OOc+00 


l.lle-16 


1 afIT 

lOE 1 ^ 


9.24e-28 


3.56e-25 


O.OOc+00 


0.00e+00 


2.80e-45 


3.40e-42 


O.OOc+00 


0.00e+00 


10F 11 


1.85e-05 


1.67c-06 


0.023 


0.00e+00 


2.29e-12 


2.88c-16 


5.96e-08 


2.07e-ll 


nA RVJ 


0.002 


0.027 


8.27c-04 


8.13e-05 


0.076 


0.043 


0.173 


0.041 


nB RVJ 


1.07c-04 


0.002 


6.54c-05 


1.12c-04 


0.006 


0.008 


0.025 


0.013 


11C RVJ 


7.88e-09 


3.01c-08 


0.00e+00 


2.29c-08 


1.78c-06 


3.98e-06 


3.40e-06 


1.41e-06 


nD RVJ 


4.41e-05 


1.22c-04 


4.17c-07 


1.20c-04 


1.55e-06 


8.37c-07 


9.18e-06 


1.98c-08 


UE RVJ 


1.87c-25 


4.71e-26 


O.OOc+00 


O.OOc+00 


5.53c-33 


7.41c-34 


O.OOc+00 


O.OOc+00 


UpRVJ 


4.70c-04 


3.99e-05 


0.008 


6.20c-ll 


3.64e-17 


1.81e-19 


O.OOc+00 


2.14c-13 


11A IT 


0.040 


0.084 


0.066 


0.277 


0.003 


5.92e-04 


0.003 


0.002 


11B IT 


0.041 


0.105 


0.024 


0.146 


0.032 


0.015 


0.083 


0.044 


11C IT 


0.186 


0.003 


0.421 


0.016 


0.001 


0.003 


0.001 


0.003 


11D IT 


1.59c-05 


3.57c-05 


5.81c-04 


9.44e-04 


6.92c-10 


1.67c-10 


0.00e+00 


l.lOc-08 


11E IT 


1.65e-28 


1.02e-25 


0.00e+00 


O.OOc+00 


8.41e-45 


1.17c-41 


O.OOc+00 


O.OOc+00 


11F IT 


8.04c-07 


2.11c-08 


0.004 


9.99c-16 


2.30c-08 


2.96c-10 


1.55e-06 


2.65e-06 


12A RVJ 


5.33c-12 


3.02e-ll 


0.00e+00 


2.02e-08 


8.52e-10 


1.15c-09 


2.03e-06 


6.45c-10 


12B RVJ 


2.91c-13 


2.24c-12 


O.OOc+00 


2.87c-09 


1.44e-10 


6.93c-10 


4.17c-07 


l.OOc-10 


12C RV.I 


3.81c-18 


4.79e-18 


O.OOc+00 


6.66e-16 


2.42e-16 


4.73c-16 


0.00e+00 


4.66e-15 
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Table 6 — Continued 



Model 




NGC 3379 






NGC 4278 






^KS 


1 Kuipcr 


P o 


Prs 


Pks 


1 Kuipcr 


P o 


Prs 


12D" ' 


1.83e-20 


7.78c-21 


0.00e+00 


0.00e+00 


7.05e-20 


8.39c-22 


0.00e+00 


0.00e+00 


12E J 


1.20e-25 


2.37e-26 


O.OOc+00 


O.OOc+00 


3.14e-31 


2.07c-29 


O.OOc+00 


O.OOc+00 


12F 


4.66c-07 


7.47c-07 


1.35e-04 


6.95e-14 


2.17e-20 


1.14c-23 


O.OOc+00 


0.00e+00 


i o A TT 

12A 11 


0.045 


0.105 


0.101 


0.116 


3.35e-04 


7.88e-05 


0.005 


2.76e-04 


12B 


0.454 


0.653 


0.148 


0.770 


0.150 


0.007 


0.121 


9.39c-04 


12C 


0.003 


5.98c-04 


0.140 


6.62c-04 


5.07e-09 


5.77c-09 


1.49c-06 


3.52e-08 


12D 


9.56e-07 


1.67e-06 


2.62c-04 


2.29e-05 


1.85e-13 


6.09e- 14 


O.OOc+00 


3.10e-12 


12E 


2.78e-28 


1.57o-25 


O.OOc+00 


0.00e+00 


0.00e+00 


1.14c-42 


O.OOc+00 


0.00e+00 


1 OT7IT 


8.30e-07 


9.82c-08 


0.004 


O.OOc+00 


1.39e-09 


1.71c-12 


2.15c-06 


5.84e-08 


1 o a RV T 

13A n J 


1.24e-04 


3.74c-04 


7.24e-05 


0.004 


0.001 


0.002 


0.004 


0.007 


1 ooRVJ 
13B ' 


4.02e-05 


1.31e-04 


2.20e-05 


0.003 


1.13c-05 


2.40c-05 


5.53e-05 


2.23c-04 




6.62e-13 


5.54e-13 


O.OOc+00 


4.30e-ll 


6.29c-14 


1.05c-14 


O.OOc+00 


7.27e-14 




5.35e-09 


5.10o-09 


O.OOc+00 


2.27c-ll 


3.74e-06 


4.37e-07 


4.49e-04 


1.18e-05 


1 OTlRV T 

13E 


2.68c-23 


5.91o-24 


O.OOc+00 


0.00e+00 


6.99e-35 


3.78c-33 


O.OOc+00 


0.00e+00 


13F^ V J 


1.99c-08 


2.00e-08 


2.03c-06 


5.75e-12 


4.10e-24 


1.80e-28 


0.00e+00 


0.00e+00 


1 o A TT 

13 A 


0.375 


0.208 


0.487 


0.056 


0.307 


0.196 


0.367 


0.078 




0.728 


0.903 


0.222 


0.707 


0.552 


0.477 


0.215 


0.438 


13C 


0.084 


0.099 


0.205 


0.004 


1.07e-04 


1.48&-04 


9.03c-04 


2.91e-04 


13D 


0.005 


0.010 


0.026 


0.013 


2.54e-07 


1.07c-07 


2.38c-07 


1.07e-06 


13E 1 1 


3.80e-26 


1.92c-23 


0.00e+00 


O.OOc+00 


1.06c-43 


1.38c-40 


O.OOc+00 


O.OOc+00 


13F 


2.58e-05 


7.92c-06 


0.008 


5.55e-16 


2.83c-10 


7.03c-09 


0.00e+00 


6.72e-09 


1 A A RVJ 


2.77e-07 


6.89e-07 


0.00e+00 


4.96c-10 


3.34c-04 


1.30e-04 


0.003 


6.39e-04 


1 /< t> RV J 
14B v J 


1.34e-10 


3.33c-10 


0.00e+00 


2.93e-10 


2.60c-08 


3.35c-08 


5.13c-06 


4.68c-07 


14C 


4.94e-19 


1.83c-19 


0.00e+00 


O.OOc+00 


3.61e-18 


2.11c-19 


0.00e+00 


O.OOc+00 


14D" v J 


6.87e-10 


5.52e-10 


O.OOc+00 


0.00e+00 


6.96e-06 


1.25c-06 


4.06c-04 


1.03e-10 


14E 


6.72e-27 


6.11c-28 


0.00e+00 


0.00e+00 


5.05e-37 


2.14c-35 


O.OOc+00 


0.00e+00 




4.38c-10 


3.35e-10 


2.38c-07 


0.00e+00 


3.17e-24 


3.97c-29 


O.OOc+00 


0.00e+00 


1 j a TT 
14A 1A 


0.804 


0.553 


0.351 


0.528 


0.661 


0.654 


0.302 


0.438 


14B 


0.780 


0.785 


0.291 


0.536 


0.890 


0.779 


0.462 


0.449 


14C IT 


0.049 


0.045 


0.171 


0.018 


3.24e-05 


4.32e-05 


3.29e-04 


3.36e-04 


14D IT 


1.70c-04 


3.22e-04 


0.001 


0.002 


1.14c-08 


5.67c-09 


0.00e+00 


2.76c-08 


14E IT 


3.82c-28 


3.38c-25 


O.OOc+00 


O.OOc+00 


O.OOc+00 


O.OOc+00 


0.00e+00 


0.00e+00 


14F IT 


1.95e-06 


5.16c-06 


0.001 


O.OOc+00 


2.21c-08 


3.94c-09 


1.31c-06 


4.57c-07 


15A RVJ 


4.20c-06 


8.93c-06 


1.79c-07 


1.53c-09 


0.014 


0.003 


0.012 


0.002 


15B RVJ 


2.07c-09 


4.54c-09 


O.OOc+00 


1.66&-09 


1.35e-07 


1.48e-07 


1.52c-05 


5.53e-07 


15C RVJ 


8.30e-18 


3.25c-18 


O.OOc+00 


O.OOc+00 


1.61c-15 


1.29c-16 


0.00e+00 


O.OOc+00 


15D RVJ 


2.31e-ll 


1.46e-ll 


O.OOc+00 


O.OOc+00 


1.32e-05 


3.34c-06 


3.20c-04 


8.85c-07 


15E RVJ 


3.68e-25 


4.02c-26 


O.OOc+00 


O.OOc+00 


5.70c-35 


1.48e-36 


O.OOc+00 


0.00e+00 


15F RVJ 


2.64c-10 


1.92c-10 


1.19c-07 


1.74c-14 


3.81e-24 


4.03c-29 


0.00e+00 


0.00e+00 


15A IT 


0.193 


0.303 


0.199 


0.065 


0.001 


0.004 


0.002 


2.36e-04 


15B IT 


0.871 


0.647 


0.366 


0.444 


0.837 


0.509 


0.483 


0.066 


15C IT 


0.298 


0.151 


0.454 


0.035 


2.92c-04 


5.57c-04 


9.97c-04 


5.96c-05 


15D IT 


0.052 


0.042 


0.194 


0.010 


2.04c-07 


2.11c-07 


2.86c-06 


2.87c-07 


15E IT 


1.72e-26 


1.51o-23 


O.OOc+00 


O.OOc+00 


1.96e-44 


2.61c-41 


O.OOc+00 


O.OOc+00 


15F IT 


4.27c-06 


1.28e-06 


0.005 


l.lle-16 


4.95e-09 


2.47e-08 


1.19c-07 


1.67c-07 



-31 - 



Table 6 — Continued 



Model 




NGC 3379 






NGC 4278 






Pks 


^Kuipcr ^~X 2 


Prs 


Pks 


-^Kuipcr ^~X 2 


Prs 



16A RVJ 
16B RVJ 
16C RVJ 
16D RVJ 

16E RVJ 

16F RVJ 

16A IT 
16B IT 
16C IT 
16D IT 
16E IT 
16F IT 

17A RVJ 
17B RVJ 
17C RVJ 
17D RVJ 
17E RVJ 
17F RVJ 

17A IT 
17B IT 
17C IT 
17D IT 
17E IT 
17F IT 

18A RVJ 
lgB RVJ 
18C RVJ 
18D RVJ 
1 8E RVJ 
18F RVJ 

18A IT 
18B IT 
18C IT 
18D IT 
18E IT 
18F IT 

lgA RVJ 
19B RVJ 
lgC RVJ 
lgD RVJ 
1 9E RVJ 
lgF RVJ 

19A IT 

19B IT 
19C IT 



1.60c-06 
1.51c-09 
3.27c-15 
1.51c-ll 
9.27c-24 
6.85e-12 
0.832 
0.585 
0.009 
1.91c-04 
8.10e-26 
2.78c-05 
3.93e-07 
3.17c-ll 
2.32e-15 
1.90e-06 
9.33e-27 
1.72c-17 
0.011 
0.201 
0.158 
7.79e-04 
3.21e-24 

0.003 
8.79c-07 
1.91c-09 
6.48e-16 
4.66c-07 
1.23e-24 
1.09e-12 
0.698 
0.904 
0.234 
0.013 
4.10e-25 
2.56c-05 
2.89e-06 
2.33c-09 
2.89c-16 
2.21c-08 
2.16e-24 
1.35e-14 
0.332 
0.348 
0.073 



2.66c-06 
2.16e-09 
1.43e-15 
8.61c-12 
1.10c- 24 
3.83c-12 
0.778 
0.629 
0.018 
3.26e-04 
5.55e-23 
9.86c-06 
3.64e-05 
4.47e-ll 
9.22c-16 
2.09e-06 
9.92e-24 
8.76c-15 
0.046 
0.225 
0.222 
0.001 
1.66e-21 
3.74e-04 
1.27c-06 
2.46e-09 
2.32c-16 
4.69c-07 
1.21e-25 
5.07c-13 
0.544 
0.771 
0.141 
0.026 
1.96e-22 
2.57c-05 
5.02e-06 
3.67c-09 
1.08c-16 
1.90c-08 
2.36e-25 
5.16e-15 
0.508 
0.546 
0.076 



3.58e-07 
O.OOc+OO 
O.OOc+00 
O.OOc+OO 
O.OOc+00 
O.OOc+OO 
0.294 
0.204 
0.044 
0.001 
O.OOc+00 

0.008 
O.OOc+OO 
O.OOc+00 
0.00e+00 
3.58c-07 
O.OOc+OO 
O.OOc+00 
0.011 
0.219 
0.174 
0.002 
0.00e+00 

0.042 
3.16c-06 
5.96c-08 
O.OOc+OO 
0.00e+00 
O.OOc+OO 
O.OOc+00 
0.254 
0.273 
0.325 
0.020 
0.00e+00 

0.004 
1.19c-07 
O.OOc+OO 
O.OOc+00 
O.OOc+OO 
O.OOc+00 
O.OOc+OO 
0.117 
0.163 
0.217 



3.50c-06 
6.68c-08 
2.22c-16 
2. lie- 15 
O.OOc+OO 
4.45e-13 
0.397 
0.471 
0.004 
0.003 
O.OOc+OO 
4.44c-16 
4.18c-ll 
2.43e-12 
O.OOc+OO 
1.27c-12 
O.OOc+OO 
0.00e+00 
0.051 
0.110 
0.073 
0.055 
O.OOc+OO 
3.73e-09 
5.40c-04 
1.49c-06 
O.OOc+OO 
6.62e-12 
O.OOc+OO 
O.OOc+OO 
0.323 
0.485 
0.026 
0.067 
O.OOc+OO 
1.52e-14 
l.OOc-ll 
1.84c-ll 
O.OOc+OO 
5.55e-16 
O.OOc+OO 
O.OOc+00 
0.538 
0.830 
0.031 



1.09c-05 
1.58e-09 
2.58c-16 
6.98e-07 
1.21c-32 
2.56e-23 
0.498 
0.276 
5.22e-07 
1.54e-09 
6.45e-44 
1.43c-10 
6.94c-05 
4.57c-ll 
2.88c-14 

0.003 
6.76c-39 
3.98c-31 
0.800 
0.192 
0.008 
0.001 
7.56c-41 
9.88e-10 
1.71c-06 
1.35c-12 
1.34c-16 
6.30e-04 
1.56e-36 
1.35e-26 
0.341 
0.739 
6.34e-04 
1.44e-05 
8.41e-44 
2.79c-08 
6.84e-04 
4.19c-09 
1.30c-15 
1.63c-04 
1.37c-35 
3.44e-25 
0.759 
0.176 
7.73e-05 



4.38c-06 
1.52e-09 
1.31c-17 
6.30c-08 
9.83c-30 
2.34e-28 
0.121 
0.129 
3.61c-07 
1.16c-10 
3.10c-41 
1.57c-13 

0.001 
2.13c-ll 
1.79c-15 
3.78c-05 
5.14c-36 
4.67c-29 
0.554 
0.561 
0.010 
9.68c-04 
1.82c-38 
9.38c-09 
2.96e-07 
3.56c-13 
5.15c-18 
l.llc-04 
2.92e-35 
3.11c-32 
0.118 
0.475 
7.53e-04 
1.06c-05 
6.59e-41 
1.63c-ll 
4.20c-04 
2.39e-09 
6.86c-17 
7.67e-06 
2.83c-37 
1.45e-30 
0.344 
0.087 
6.35e-05 



1.87e-04 
3.58c-07 
0.00e+00 
2.49e-05 
O.OOc+OO 
0.00e+00 
0.377 
0.328 
3.91e-05 
0.00e+00 
O.OOc+OO 
1.79c-07 
8.11c-05 
O.OOc+OO 
O.OOc+00 

0.046 
O.OOc+OO 
O.OOc+00 
0.443 
0.112 
0.030 
3.70c-04 
O.OOc+OO 
5.96e-08 
5.68e-05 
0.00e+00 
O.OOc+OO 

0.001 
0.00e+00 
O.OOc+OO 
0.458 
0.311 
0.007 
2.74c-06 
O.OOc+OO 
2.18c-05 

0.003 
1.19e-07 
O.OOc+OO 

0.001 
O.OOc+OO 
O.OOc+OO 
0.350 
0.125 
0.001 



1.75c-05 
6.88e-09 
0.00e+00 
1.67c-15 
0.00e+00 
0.00e+00 
0.022 
0.007 
6.95e-07 
8.82c-08 
O.OOc+OO 
3.03c-09 
9.88c-04 
1.40c-09 
O.OOc+OO 
1.22e-05 
O.OOc+OO 
O.OOc+00 
0.118 
0.155 
0.036 
7.76c-04 
O.OOc+OO 
1.98c-07 
3.95c-08 
2.93c-13 
0.00e+00 
5.70c-07 
O.OOc+OO 
O.OOc+OO 
0.068 
0.293 
0.002 
5.47c-05 
O.OOc+OO 
1.96c-07 
5.41c-04 
3.73c-08 
O.OOc+OO 
2.66c-09 
O.OOc+OO 
O.OOc+00 
0.211 
0.253 
5.00e-04 
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Model 




NGC 3379 






NGC 4278 






Pks 


^Kuipcr ^X 2 


Prs 


Pks 


^Kuiper ^~X 2 


Prs 



19D IT 


0.002 


0.004 


19E IT 


7.61c-25 


3.26e-22 


19F IT 


8.29c-05 


8.49e-05 


20A RVJ 


2.87c-08 


3.30c-08 


20B RVJ 


3.35e-ll 


3.10e-ll 


20C HVJ 


4.68e-18 


1.17c-18 


20D RVJ 


2.72c-08 


2.21e-08 


20E RVJ 


2.23e-25 


1.89c-26 


20F RVJ 


2.42e-15 


7.59e-16 


20A IT 


0.648 


0.451 


20B IT 


0.902 


0.934 


20C IT 


0.217 


0.268 


20D IT 


0.026 


0.048 


20E IT 


4.39c-23 


1.93e-20 


20F IT 


5.12e-05 


2.85e-05 


21A RVJ 


1.90e-ll 


1.96e-ll 


21B RVJ 


4.57c-13 


5.07e-13 


21C RV.J 


2.23c-20 


5.38e-21 


21D RVJ 


1.77e-16 


5.25e-17 


21E RVJ 


3.49c-26 


3.16e-27 


21F RVJ 


3.13c-14 


1.30e-14 


21A IT 


0.017 


0.124 


21B IT 


0.021 


0.046 


21C IT 


0.013 


0.017 


21D IT 


5.24e-04 


9.19e-04 


21E IT 


1.14e-27 


9.70c-25 


21F IT 


0.014 


1.88c-04 


22A RVJ 


8.20e-12 


3.49e-10 


22B RVJ 


7.93e-13 


9.61e-13 


22C RVJ 


1.03c-19 


2.40e-20 


22D RVJ 


9. 69c- 14 


3.88c-14 


22E RVJ 


6.13c-29 


2.61e-27 


22F RVJ 


7.98c-14 


1.95c-12 


22A IT 


0.184 


0.566 


22B IT 


0.887 


0.921 


22C IT 


0.424 


0.122 


22D IT 


0.041 


0.020 


22E IT 


6.17c-25 


2.35e-22 


22F IT 


2.66e-04 


6.07c-05 


23A RVJ 


1.95c-13 


1.66e-13 


23B RVJ 


5.78c-15 


5.19c-15 


23C RVJ 


1.23c-22 


2.17e-23 


23D RVJ 


2.48c-18 


5.66e-19 


23E RVJ 


8.20e-27 


6.84e-28 


23F rvj 


7.02e-17 


2.02e-17 



0.011 


0.129 


6.59e-06 


O.OOc+00 


O.OOc+00 


1.50e-42 


0.006 


7.14c-12 


6.05c-09 


5.96c-08 


1.60c-07 


6.44c-07 


O.OOc+00 


8.23c-10 


1.76c-13 


O.OOc+00 


0.00e+00 


1.96e-19 


O.OOc+00 


O.OOc+00 


3.82c-04 


O.OOc+00 


0.00e+00 


1.32e-35 


O.OOc+00 


0.00e+00 


3.49c-27 


0.294 


0.302 


0.038 


0.326 


0.568 


0.320 


0.244 


0.039 


0.008 


0.041 


0.271 


0.002 


O.OOc+00 


O.OOc+00 


1.32e-40 


0.005 


2.00c-12 


2.72c-06 


0.00e+00 


0.00e+00 


3.08c-09 


O.OOc+00 


3.16c-13 


2.07c-12 


0.00e+00 


O.OOc+00 


7.10c-21 


0.00e+00 


O.OOc+00 


9.68c-13 


O.OOc+00 


0.00e+00 


7.65e-33 


O.OOc+00 


0.00e+00 


9.56e-26 


0.003 


0.331 


0.021 


0.010 


0.051 


0.356 


0.027 


0.003 


2.22e-04 


0.002 


0.049 


1.36e-06 


O.OOc+00 


O.OOc+00 


3.05e-43 


0.125 


2.44c-10 


2.41e-13 


0.00e+00 


3.84e-ll 


5.64e-ll 


0.00e+00 


9.41e-10 


1.80e-14 


O.OOc+00 


O.OOc+00 


3.25e-22 


O.OOc+00 


O.OOc+00 


5.84e-08 


0.00e+00 


O.OOc+00 


8.32e-43 


0.00e+00 


O.OOc+00 


2.91c-33 


0.095 


0.606 


0.045 


0.331 


0.793 


0.845 


0.408 


0.046 


0.003 


0.160 


0.090 


3.17c-05 


0.00e+00 


O.OOc+00 


3.48c-43 


0.005 


2.66e-ll 


1.46e-09 


O.OOc+00 


0.00e+00 


2.18c-09 


O.OOc+00 


l.lle-16 


4.42e-13 


0.00e+00 


0.00e+00 


1.08c-22 


0.00e+00 


0.00e+00 


1.22e-12 


O.OOc+00 


0.00e+00 


7.58e-33 


0.00e+00 


0.00e+00 


1.87e-26 



2.70e-07 


2.38e-06 


1.72e-05 


7.00e-40 


O.OOc+00 


O.OOc+00 


1.02c-10 


4.77c-07 


3.89c-07 


1.70e-07 


3.26c-05 


1.71c-06 


5.14c-14 


O.OOc+00 


1.69c-12 


1.00e-20 


0.00e+00 


O.OOc+00 


2.76c-05 


2.14e-04 


3.98c-08 


2.78c-34 


0.00e+00 


O.OOc+00 


4.88e-33 


O.OOc+00 


O.OOc+00 


0.005 


0.221 


0.005 


0.131 


0.434 


0.044 


0.009 


0.050 


0.021 


8.21e-04 


0.002 


0.003 


1.25e-37 


O.OOc+00 


O.OOc+00 


4.21c-09 


1.93c-04 


8.22c-05 


1.96c-10 


2.84e-05 


4.71e-09 


6.71c-13 


2.98e-07 


3.46e-12 


1.62e-22 


0.00e+00 


O.OOc+00 


5.74c-14 


O.OOc+00 


O.OOc+00 


3.24e-32 


O.OOc+00 


O.OOc+00 


4.10c-31 


0.00e+00 


O.OOc+00 


0.145 


0.004 


0.158 


0.135 


0.477 


0.059 


2.03e-04 


0.003 


2.71c-04 


2.17c-07 


4.17c-07 


5.51c-06 


2.36e-40 


O.OOc+00 


O.OOc+00 


1.83c-ll 


O.OOc+00 


3.04c-08 


6.61c-09 


O.OOc+00 


1.69c-08 


2.44c-14 


O.OOc+00 


1.19e-13 


5.16c-25 


O.OOc+00 


O.OOc+00 


3.93c-12 


1.85e-04 


O.OOc+00 


1.06e-39 


O.OOc+00 


O.OOc+00 


5.48e-34 


O.OOc+00 


O.OOc+00 


0.103 


0.035 


0.358 


0.575 


0.442 


0.426 


7.28e-04 


0.025 


0.019 


9.40e-06 


2.01e-04 


6.87c-04 


2.40e-40 


O.OOc+00 


O.OOc+00 


1.13c-08 


O.OOc+00 


4.51e-06 


3.55c-10 


2.38e-06 


3.86c-09 


1.19c-13 


0.00e+00 


1.34e-12 


1.64e-24 


O.OOc+00 


O.OOc+00 


3.14c-15 


O.OOc+00 


O.OOc+00 


2.16c-34 


0.00e+00 


O.OOc+00 


5.39e-32 


0.00e+00 


0.00e+00 
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NGC 3379 



NGC 4278 



Kuipcr 



Prs 



Pks 



Kuipcr 



-Prs 



Model 

Pks 



OQ A IT 
ISA. 


0.011 


0.022 


zoB 


0.191 


0.344 


23C 


0.004 


0.002 


oqt^IT 

zoJJ 


1.37e-05 


o no,-, nc; 

z.Uoe-Uo 




3.08e-25 


2.64e-22 


oq-pIT 


o a n a 
z.4Ue-U4 


1.38e-04 


o/i A RVJ 
Z4A 


2.57e-09 


o.bze-uy 


Z4r> 


y.yoe-lZ 


1 i n« 1 1 


Z4L- 


i.uye-iy 


o c^ft^, on 
z.ooe-zU 


Z4U 


Q AOn 1 Q 

o.UZe- lo 


l.zye-io 


O/irpRVJ 
Z4rA/ 


z. / /e-zY 


1.95e-28 


o /it? RVJ 
Z4r 


1.65e-12 


7.90c-13 


0/1 A IT 
Z4A 


n nftQ 


n riQn 


Z4r> 


0.312 


0.495 


z4L/ 


0.060 


0.110 


o/ir~\IT 


4.8UC-U4 


7.97e-U4 


z4r!i 


cr o^n O/l 


2.25e-21 


O/ITpIT 

24b 


1.60c-04 


1.53e-04 


or, A RVJ 
zoA 


4. 19e-15 


1.08e-12 


OCR RVJ 

ZorS 


O Ol 1 ft 
Z.ZlC-lD 


1 Qflr. 1 E 

i.oUe-io 


ZOL/ 


i fti on 


1 /I on 


ot-riRVJ 
zoJJ 


1.90e-12 


8.34e-13 


or-cRVJ 

ZDrir 


2.59e-27 


2.52e-24 


ok TP RVJ 
ZOr 


6.34e-22 


4.26c-19 


orr A IT 
ZoA 


q ro„ nft 

o.ooe-Uo 


1 1 O^ H/1 


okt^IT 
ZoB 


3.30e-05 


c o n „ n cr 

D.oUc-Uo 


z5C 


5.03e-ll 


3.89e-ll 


25D IT 


3.45c-13 


1.71c-13 


25E IT 


1.68c-20 


8.88e-18 


25F IT 


8.84c-09 


1.58c-07 


26A RVJ 


7.05c-16 


2.51e-13 


26B RVJ 


3.61c-15 


2.97e-15 


26C RVJ 


1.20c-20 


2.41e-21 


26D RVJ 


3.74c-13 


1.61e-13 


26E RVJ 


1.19c-27 


6.74e-25 


26F RVJ 


1.92c-16 


5.57c-15 


26A IT 


9.89e-05 


3.96e-04 


26B IT 


0.001 


0.003 


26C IT 


3.73c-08 


4.34e-08 


26D IT 


1.69c-10 


1.27e-10 


26E IT 


2.13c-19 


6.54e-17 


26F IT 


0.001 


0.005 


27A RVJ 


4.11c-15 


1.87c-15 


27B RVJ 


1.34c-17 


4.94c-18 


27C RVJ 


1.38e-23 


1.58e-24 
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0.003 


0.042 


0.149 


0.049 


0.399 


0.281 


0.037 


0.013 


1.23c-06 


2.10e-04 


0.001 


2.74c-08 


O.OOc+00 


0.00e+00 


6.81c-41 


0.008 


1.05c-ll 


2.76c-10 


5.96c-08 


3.83c-08 


2.20e-06 


O.OOc+00 


6.32c-09 


1.55e-13 


O.OOc+00 


O.OOc+00 


1.56c-20 


O.OOc+00 


0.00e+00 


2.46e-07 


0.00e+00 


O.OOc+00 


5.21c-41 


0.00e+00 


0.00e+00 


4.65e-31 


0.027 


0.183 


0.189 


0.090 


0.516 


0.270 


0.063 


0.113 


5.66c-04 


0.004 


0.037 


9.62c-06 


0.00e+00 


O.OOc+00 


6.40e-41 


0.006 


2.87c-09 


2.03c-07 


0.00e+00 


2.27c-13 


2.91c-15 


O.OOc+00 


6.66c-16 


3.75c-17 


O.OOc+00 


O.OOc+00 


9.06c-20 


0.00e+00 


O.OOc+00 


1.64e-08 


0.00e+00 


O.OOc+00 


1.20e-35 


0.00e+00 


O.OOc+00 


1.20e-29 


3.58e-07 


1.10e-10 


0.019 


4.83c-06 


1.51c-05 


4.67c-04 


O.OOc+00 


1.34e-09 


2.49c-12 


0.00e+00 


l.lle-16 


4.20e-10 


O.OOc+00 


O.OOc+00 


5.39e-28 


1.73c-06 


4.08c-06 


1.74e-18 


0.00e+00 


O.OOc+00 


4.06c-12 


O.OOc+00 


4.44e-16 


8.48c-15 


O.OOc+00 


O.OOc+00 


2.44c-19 


O.OOc+00 


O.OOc+00 


8.89c-08 


O.OOc+00 


0.00e+00 


2.51e-37 


0.00e+00 


O.OOc+00 


3.38c-29 


1.07e-05 


3.15e-05 


0.004 


6.83c-04 


0.002 


0.001 


2.98c-07 


4.49e-06 


5.10c-ll 


O.OOc+00 


4.09e-09 


2.93c-10 


O.OOc+00 


O.OOc+00 


4.74e-29 


0.033 


2.23c-05 


1.19e-13 


O.OOc+00 


O.OOc+00 


2.40e-14 


O.OOc+00 


0.00e+00 


5.62e-18 


0.00e+00 


O.OOe+00 


7.95e-26 



0.026 


0.123 


0.008 


0.060 


0.213 


0.010 


7.79c-07 


3.41c-05 


1.32c-05 


2.50e-09 


2.98c-07 


4.97c-07 


4.79c-38 


0.00e+00 


O.OOc+00 


1.72c-12 


5.96e-08 


2.71c-08 


5.51c-07 


3.23e-05 


4.95c-06 


1.07c-13 


O.OOc+00 


1.40c-12 


l.Olc-24 


O.OOc+00 


O.OOc+00 


4.06c-12 


2.32e-04 


O.OOc+00 


4.85e-41 


0.00e+00 


O.OOc+00 


2.96e-37 


O.OOc+00 


O.OOc+00 


0.018 


0.462 


0.003 


0.129 


0.229 


0.245 


4.91c-04 


0.002 


0.007 


4.75e-06 


6.02e-06 


5.93e-06 


3.97c-38 


O.OOc+00 


O.OOc+00 


4.65c-10 


3.22e-06 


l.lOc-05 


4.19c-13 


O.OOc+00 


5.58e-13 


2.92c-16 


O.OOc+00 


O.OOc+00 


1.42e-21 


O.OOc+00 


O.OOc+00 


1.32c-10 


2.98e-07 


2.60c-09 


1.20e-32 


O.OOc+00 


O.OOc+00 


8.92e-27 


O.OOc+00 


O.OOc+00 


0.086 


0.017 


0.004 


1.27c-05 


0.015 


7.42c-06 


3.01c-13 


O.OOc+00 


1.32c-ll 


4.44c-13 


O.OOc+00 


5.79c-12 


2.42e-25 


0.00e+00 


O.OOc+00 


3.80c-16 


0.00e+00 


O.OOc+00 


5.98c-10 


0.00e+00 


2.83c-14 


4.45c-14 


O.OOc+00 


6.66c-16 


5.50c-21 


O.OOc+00 


O.OOc+00 


1.74c-09 


5.96c-07 


3.10c-08 


2.59e-34 


O.OOc+00 


O.OOc+00 


2.50e-28 


O.OOc+00 


O.OOc+00 


0.003 


0.013 


0.002 


1.48e-06 


0.048 


3.62e-06 


9.79c-12 


0.00e+00 


6.68c-10 


1.91c-ll 


0.00e+00 


1.20c-10 


2.26e-26 


O.OOc+00 


O.OOc+00 


2.64e-13 


O.OOc+00 


8.80e-13 


4.67c-16 


0.00e+00 


2.22c-16 


1.36c-18 


O.OOc+00 


O.OOc+00 


3.92e-28 


0.00e+00 


0.00e+00 
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Table 6 — Continued 



Model 




NGC 3379 






NGC 4278 






Pks 


^~*Kuipcr 






Pks 


-^Kuipcr 






27D RV.I 


2.69c-19 


4.60e-20 


0.00e+00 


O.OOc+00 


9.91c-16 


2.20c-19 


0.00e+00 


0.00e+00 


27E RVJ 


4.33e-27 


2.88c-28 


0.00e+00 


O.OOc+00 


1.24c-31 


6.52e-31 


0.00e+00 


0.00e+00 


27F RVJ 


4.72e-22 


5.95e-23 


0.00e+00 


O.OOc+00 


4.99e-29 


2.99e-35 


0.00e+00 


O.OOc+00 


27A IT 


6.35e-06 


8.72c-06 


2.15e-06 


5.66c-10 


0.083 


0.006 


0.084 


0.006 


27B IT 


1.64e-05 


3.06c-05 


9.30c-06 


2.62e-05 


0.002 


1.75c-05 


0.014 


7.18c-05 


27C IT 


4.54e-12 


2.95e-12 


0.00e+00 


6.01e-ll 


1.76c-13 


1.55e-14 


O.OOc+00 


7.22e-14 


27D IT 


2.66e-15 


9.48e-16 


0.00e+00 


O.OOc+00 


5.39c-15 


1.42c-17 


0.00e+00 


O.OOc+00 


27E IT 


5.43e-20 


2.00e-20 


O.OOc+00 


O.OOc+00 


8.08e-27 


1.34e-25 


0.00e+00 


O.OOc+00 


27F IT 


1.08c-07 


1.06c-07 


1.36c-05 


2.04e-05 


4.21c-16 


4.71c-20 


O.OOc+00 


O.OOc+00 


28A RVJ 


3.81c-13 


3.68e-13 


O.OOc+00 


O.OOc+00 


5.83c-09 


7.29c-10 


4.77c-07 


8.57e-13 


28B RVJ 


5.01e-15 


4.98e-15 


O.OOc+00 


O.OOc+00 


1.14c-14 


3.07c-15 


0.00e+00 


9.66e-15 


28C RVJ 


2.61c-20 


6.54e-21 


0.00e+00 


O.OOc+00 


1.78c-19 


5.48e-21 


0.00e+00 


O.OOc+00 


28D RVJ 


3.94e-14 


1.62e-14 


0.00e+00 


O.OOc+00 


4.29e-09 


1.82c-ll 


0.00e+00 


7.06c-10 


28E rvj 


6.94e-27 


5.83e-28 


O.OOc+00 


O.OOc+00 


1.04c-34 


1.14c-36 


O.OOc+00 


0.00e+00 


28F RVJ 


2.91e-14 


1.21c-14 


O.OOc+00 


O.OOc+00 


5.29c-27 


1.66c-32 


O.OOc+00 


0.00e+00 


28A IT 


1.08c-05 


1.65e-05 


5.84e-06 


7.18c-05 


0.003 


4.00c-04 


0.007 


8.08c-04 


28B IT 


0.001 


0.003 


4.16e-04 


0.002 


0.003 


1.24e-04 


0.028 


4.04e-04 


28C IT 


6.38e-10 


6.21c-10 


0.00e+00 


3.07c-08 


2.45e-13 


3.26c-14 


O.OOc+00 


9.36e-13 


28D IT 


9.54e-14 


4.64e-14 


O.OOc+00 


1.01e-13 


1.84e-14 


1.83c-16 


O.OOc+00 


O.OOc+00 


28E IT 


6.92e-21 


2.39e-21 


O.OOc+00 


O.OOc+00 


3.34e-32 


2.51c-31 


O.OOc+00 


0.00e+00 


28F IT 


5.66c-04 


8.88c-04 


0.023 


2.58c-06 


3.18c-16 


7.75c-21 


O.OOc+00 


2.22c-16 
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Fig. 4. — Model XLF for all the different LMXB population models listed in Tables [3] 
and HI In each panel all the PS parameters are kept constant, and only the modeling of 
transient systems changes. For comparison the observed XLFs of NGC3379 (dark red) and 
of NGC4278 (orange) are drawn. 
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Fig. 5. — Analysis of the LMXB population for complete list of our models. We show the 
contribution of different sub-populations by separating the LMXBs into groups of systems 
with different donor and accretor stellar types. 
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